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I.  INTRODUCTION 


Contract  AF  33(657)-10567  provides  for  an  investigation  of  the  recently  developed 
superconducting  magnets  and  their  application  to  MHD  power  generators  for  space  use. 

The  work  will  include  theoretical  analysis  of  an  MHD  generator  with  particular  emphasis 
on  the  magnetic  field  requirements  and  a  study  of  the  problems  of  superconducting  magnet 
design  to  fulfill  these  requirements. 

To  attain  this  end,  models  of  a  superconducting  magnet,  cooling  equipment  and  an 
MHD  generator  will  be  analyzed.  The  total  power  system  including  a  reactor  heat  source 
and  radiator  will  then  be  optimized  on  a  power  to  weight  basis  to  find  the  optimum  magnetic 
field.  With  this  completed,  a  comparison  will  then  be  made  between  superconducting 
magnets  and  conventional  magnets  to  illustrate  the  advantages  or  disadvantages  of  both. 

• 

II.  MHD  POWER  GENERATOR 


A.  General  Considerations 

The  purpose  of  the  MHD  generator  model  for  this  study  is  to  provide  an  estimate  of  the 
magnetic  field  requirements  which  must  be  satisfied  by  superconducting  magnets  for  space 
power  generation.  Particular  areas  of  study  are;  heat  transfer  and  wall  friction  as  they  are 
affected  by  magnetic  fields;  insulation  requirements  to  reduce  the  heat  conducted  to  the 
superconducting  magnet  from  the  hot  plasma  or  high  temperature  walls;  effects  due  to  high 
fields  such  as  ion  slip,  losses  such  as  end  losses,  and  joule  heating;  and  the  flux  homogeneity 
required  in  the  channel  to  prevent  circulating  currents  and  other  effects. 

A  power  unit  to  fulfill  this  requirement  will  peed  to  have  an  electrical  output  of  at 
least  1  MW,  be  representative  of  the  state  of  the  art  in  materials  and  technology  and  have 
a  reasonable  operating  efficiency.  It  must  be  realized  that  MHD  units  are  competing  with  a 
much  more  advanced  state  of  the  art  in  turbine-generator  power  systems.  Therefore  sizable 
advantages  of  these  units  must  be  shown  to  make  their  development  economically  feasible. 


2 


Investigators^  have  shown  that  the  Rankine  cycle  Is  superior  to  the  Brayton  eycle 
for  space  power  generation.  This  is  largely  due  to  the  radiator  size  which  results  from  the 
requirement  of  heat  rejection  from  a  gas  at  a  relatively  low  temperature  characteristic  of 
the  Brayton  cycle.  Selecting  a  Rankine  cycle  fluid  with  known  thermodynamic  and  electrical 

properties  leads  to  many  serious  problems  and  compromises.  Rankine  cycle  fluid  requirements 

(2) 

are  discussed  by  Keenan'  '  and  possible  fluids  to  fulfill  these  requirements  are  outlined  by 
T„®. 


To  date  most  of  the  plasma  diagnostic  work  has  been  performed  on  either  combustion 
systems  or  on  gas  systems  seeded  with  an  alkali  vapor.  Electrically  the  alkali  vapors  are 
advisable  due  to  the  low  ionization  potential  and  the  knowledge  of  their  electrical  and 
thermodynamic  properties.  However/  thermodynamically  other  metal  and  salt  vapors  may 
be  more  desirable. 

Lithium  appears  to  be  the  best  fluid  because  of  its  low  molecular  weight  but  more 
important,  its  vapor  pressure  allows  reasonable  operation  in  the  temperature  range  of  1400 
to  2200°  K  which  is  the  range  of  interest  for  this  application.  It  is  still  advisable  to  use 
cesium  as  a  seeding  agent  to  obtain  reasonable  conductivities  and  maximum  power  densities. 
Power  density  will  be  maximized  later  in  this  report  for  varying  seeding  ratios.  It  will  still 
be  necessary  to  assume  some  degree  of  non- thermal  ionization  either  from  an  electron  bom¬ 
bardment  ionization  scheme  or  from  some  other  source. 

B.  MHD  Generator  Model 

Models  are  available  in  the  literature  for  various  MHD  generator  configurations.  The 

(4  5) 

one  selected  was  a  constant  Mach  number  solution  by  Swift-Hook'  '  '  since  this  analysis 
includes  both  heat-transfer  and  frictional  effects  which  are  among  the  major  losses  in  an  MHD 
generator.  The  solution  is  based  upon  maximum  electrical  power  extraction  per  unit  length 
and  describes  a  segmented  electrode  Faraday  generator  operating  at  constant  magnetic  field 
and  with  a  constant  loading  parameter  (  r/vV).  It  is  basically  applicable  to  the  lithium  vapor 
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MHD  power  generator  since  for  this  case  the  mobility  is  rather  low  and  with  reasonable 
magnetic  fields  the  Hall  parameters  p  B  will  be  in  the  range  of  1  to  20. 

Figures  1  and  2  show  the  ranges  of  applicability  of  solutions  based  on  continuous 
electrode  Faraday  and  Hall  MHD  generators  compared  to  the  segmented  electrode  Faraday 
generators. 

The  walls  of  a  duct  containing  a  condensable  vapor  near  saturation  conditions  must 
be  maintained  at  a  very  high  temperature  to  prevent  condensation  which  would  cause  short 
circuits.  This  can  be  accomplished  by  suitable  insulation  and  heating  due  to  the  stagnation 
temperature  of  the  boundary  layer,  in  general  this  temperature  will  be  100  to  200°  K  above 
the  static  temperature  of  the  gas  stream.  Neglecting  heat  transfer  and  including  frictional 
effects  causes  the  solution  to  deviate  from  the  small  perturbations  assumed  by  Swift- Hook 
but  also  leads  to  a  more  definable  solution.  It  will  be  derived  in  this  report  beyond  the  point 
where  it  deviates  from  the  solution  in  reference  5. 

The  momentum  and  energy  equations  respectively  appear  as  equations  7  and  8  in 
reference  5: 


where  ^-S  ■  is  the  frictional  force  per  unit  volume  and  ^s--  is  convective  heat  loss  per 
unit  volume  which  will  be  assumed  negligible  for  this  case. 

In  an  ionized  fluid  the  frictional  effects  are  a  function  of  the  applied  magnetic  field, 
however,  after  examining  several  of  the  solutions  giving  frictional  coefficients  as  function 
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of  magnetic  field,  it  was  determined  that  the  frictional  coefficient  assumed  by  Swift-Hook 
is  a  reasonable  assumption  and  gives  comparative  results  for  the  assumed  turbulent  flow, 
it  is: 


f 


0.046 


(3) 


where  forsc|uare  ducts  he  gives: 


Rey 


(4) 


therefore: 


N  =  1/2  f/>u2  (5) 

Now  using  the  equation  of  state,  the  Mach  number  equation  and  Ohm's  law  current 
equation: 


j  =  ^  (u  B  -  Ey) 

and  solving  equations  1  and  2  simultaneously  leads  to: 


4f- 


where 


v 


(1  -  K)  (y  -  1)  M 


0 


-1 


(6) 


(7) 


(8) 
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and 


(9) 


which  is  the  customary  loading  parameter. 


Now  integration  results  in: 


he  sets: 


giving: 


P1 


where 


and 


and  other  thermodynamic  properties  are  given  as: 


(10) 


(ID 

(12) 


(13) 


(14) 


(15) 


6 


o 

u. 


,1/2 


04) 


=  a 


,V2-fi 


(17 


It  remains  to  find  these  variables  as  a  function  of  the  independent  parameter  x  which  can 
be  found  using  the  energy  equation.  This  results  in,  after  some  degree  of  manipulation. 


d  t 

TT 


where 


x  - 
o 


=  (J£_)  (_!!_}  (_LL)  JL 
W  V  r  \ 

Pw, 


K  (1-K)  u*  B2 


Total  inlet  power  per  unit  area  is  designated  by; 
P 


w 


1 


/\  u.  C  T 
1  1  P  °i 


08) 


(19) 


(20) 


The  power  law  approximation  for  conductivity  holds  reasonably  well  for  this  case  also: 


where 


and 


ar 


-z 


y  = 


Z  - 


d  (ln<r  ) 


sTflfTTT 


p  =  constant 


d(ln^) 


d  (In  p ) 


(21) 


(22) 


(23) 


T  =  constant 


e 
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Equation  (17)  becomes  by  substituting  the  variables  as  f  (  a,  t ): 


d  t 

rr 


z  +  1  . 

_a _  t  (1  -  w) 


(24) 


where 


w  =  (z  +  1)  (3  -  y  -  0.5 


(25) 


To  this  point  the  solution  was  taken  from  Swift-Hook,  the  only  difference  being  the 
neglecting  of  the  heat  transfer  term  in  each  of  the  equations.  To  solve  this  for  the 
frictional  case,  the  coefficient  a  can  be  re-written  as  f  (  a,  t )  by  substitution  of 
the  variables  in  that  form  into  equation  11. 


,  _  f*  -0.4  .  2/3  -  1  +  n 

I na  5  /  -Cd  t  — - - - 


dx 


(26) 


where 


T-  -2n'  fi  M  (r-  D 

- TVi - 


(27) 


We  have  assumed  that  viscosity  is  represented  as  by  Swift-Hook: 


—  T11 

=  m  T 


(28) 


where 


m  =  1.09x10 


(29) 


and 

n.  =  0.44  # 

for  the  case  of  cesium  seeded  lithium  in  the  1600  to  2000°  K  temperature  range. 


(30) 
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Differentiating  equation  (26)  gives: 


d  a 


=  -Tr  a 


0.6 


20-  0.56 


dividing  (24)  by  (31)  and  separating  variables  gives: 

‘t 

/  2  -0. 56-5+5w,. 

£x  /  .  - 5 -  d  t  = 

0  ./}  t  5 


0.4 +  z  . 
a  da 


which  integrates  to: 


•  E+  tu¥- 10 ,-4+z 


where 


7r*  (i.4  +  z) 

g  _  o _ 

*  _  (p 

^  -  2^  -  0.56  +  5  w 


substituting  (33)  into  (24): 

jT  ♦  f  <*r-  if1  u*z)  t (w ' '>  d t  =  -ii 

O 

(5) 

which  is  the  form  of  the  incomplete  beta  function  '  .  This  function  is  written  asj 

B  (I  q)  =/  y(l”1)  (1  -y)(q_1)  dy 

'  Jo 

This  form  is  obtained  by  use  of  the  substitution: 


(31) 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 


(38) 
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where  y  is  the  new  variable,  equation  (36)  becomes: 


and 


X  2  d-y)  (cr,>y  <l_1>  d y - ^V'  O-rt^V (M) 

“  -  dx=-^x 

Jo 


dy 


1  +2 


<f>  =  —it - f-  g  i  —  n  -  £ ) (  1,4  +  2  ?  } 


-(-!)—■  (i  -  g ) 

O  y 


(39) 


where 


_  ,  1  +  z 

^  1.4  +  2 


(40) 


and 


w 


(41) 


In  general,  the  tables  for  the  incomplete  beta  function  are  compiled  over  increments  too 
large  for  use  and  the  values  for  1  and  q  can  occur  outside  of  the  tables.  '  Therefore,  a 
rapidly  convergent  series  seems  the  best  method  of  representation,  The  most  rapidly  con¬ 
vergent  series  representation  is: 

n  +  1 


B  (1,  q)  =  (q‘  V  (-  1)"  JL-r 

n  =  o  '  n  '  '  '  n+ 


(42) 


where  (  ^  n  ^ )  represents  the  binomial  coefficient. 


To  obtain  an  idea  of  efficiency  for  this  machine  we  will  use  the  turbine  or  isentropic  efficiency. 


f)  -  Actual  enthalpy  drop 

/fcs  Isentropic  enthalpy  drop 


(43) 
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which  result*  in: 


n*  - 


i  - 1 


i  - 


t1/ 


(44) 


this  leaves  the  only  undefined  parameter  as  the  Hall  voltage: 


E  = 


u  B2(  K  - 


x  ».  ■»  V  *  -  1  ) 


(45) 


Power  output  as  a  function  of  temperature  ra^>  for  the  case  of  negligible  heat  transfer 


is: 


d  W 

e 

dx 


=  I 


E  A  =  f  u  A  -4-  (C  T  +  -£-) 

d  x  '  p  n 


(46) 


integrating: 


W  =  Q.^ 


(47) 


where  Qj  is  the  total  inlet  convective  enthalpy: 


Q,  =  t 


,  u,  A.  C  T 

1  1  1  p  o. 


P  Ai 
W1  1 


(48) 


and 


%  =  1-t 


(49) 


Power  output  can  now  be  determined  as  a  function  of  length  by  equation  (39). 
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C.  Gos  Properties 

With  the  duct  parameters  determined  it  remains  to  define  the  gas  properties  and  the 
operating  conditions  of  the  generator.  To  obtain  the  maximum  power  density  for  the  machine 
it  is  necessary  to  determine  an  optimum  seeding  ratio.  For  a  segmented  Faraday  generator  the 
power  density  per  unit  volume  can  be  written,  over  the  range  of  interest,  as: 

=  <^u2  B2  K  (1  -  K)  (50) 

Now  assuming  the  gas  will  have  only  a  smalT  degree  of  ionization  conductivity  can 
be  written  for  the  electron-neutral  collision  processes  as: 

<?U=  ne  e  (51) 

Electron  collision  cross  section  data  for  lithium  vapors  is  obtained  from  reference  6 
and  it  is  assumed  that  for  reasonable  cesium  seeding  rates  the  mobility  is  unaffected  by  the 
addition  of  cesium.  Figure  3  shows  electron  mobility  as  a  function  of  gas  temperature. 

The  two  parameters  of  equation  (50)  which  will  determine  power  density  and  are  a 
function  of  seeding  ratio  are  conductivity  and  velocity.  Conductivity  will  increase  with 
increased  cesium  concentration  and  velocity  will  decrease  due  to  the  increasing  molecular 
weight.  Therefore  equation  (50)  will  reach  a  maximum.  To  determine  this  maximum  we 
can  write  each  of  these  variables  as  a  function  of  seeding  ratio:  ® 


2  _  i  +  e 

u  ~  “rarnccr  t 

1 1  cs 


(52) 


and 


(53) 
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by  virtue  of  the  changing  electron  density  bated  on  the  Saha  equation.  We  can  tee  that 
this  expression  is  a  weak  function  of  temperature. 


Substituting  these  expressions  into  (46)  results  in: 

:’/2 


1  +  e 


cs 


6' 

I  e 
exp  cs 

2TT 


777  e 

exp  ti 

TFT 


(54) 


This  expression  is  represented  by  figure  4  for  the  minimum  operating  temperature  of 
the  generator  and  results  in  a  rather  fiat  curve  with  a  maximum  near  0.07  which  will  be  used 
for  future  calculations. 


Knowing  the  gas  mixture,  the  average  molecblar  weight,  specific  heat,  viscosity,  and 
conductivity  can  now  be  calculated.  The  data  for  these  properties  were  obtained  from 
reference  8.  Equilibrium  values  were  assumed  near  the  upper  operating  temperature  since  the 
gas  will  have  reached  an  equilibrium  passing  through  the  reactor  heat  source  and  will  be  semi- 
frozen  on  passing  at  high  velocity  through  the  generator.  Conductivity  is  given  as. figures  5 
and  6  and  viscosity  as  figure  7. 


D.  Operating  Conditions 

When  operating  with  a  condensing  vapor  MHD  power  generation  cycle,  numerous 
considerations  must  be  made.  First  and  foremost,  it  is  advisable  to  operate  as  near  saturation 
conditions  as  possible  to  obtain  the  maximum  Carnot  efficiency.  The  upper  temperature  limit 
is  set  by  the  material  considerations  of  the  reactor  and  a  minimum  temperature  is  set  by  the 
radiating  temperature  of  the  radiator.  The  Carnot  efficiency  is  reduced  by  a  factor  of  0.58 
if  the  maximum  temperature  is  reduced  only  from  2000°  K  to  1900°  K  as  is  illustrated  by 
figure  8. 

Choosing  the  inlet  conditions  to  be  very  near  saturated  vapor  at  2000°  K  gives  a 
pressure  of  9. 2  atm  and  the  conductivity  becomes  very  low  (0.45  mho/m).  It  becomes  very 
obvious  that  for  a  lithium  vapor  cycle  some  means  of  non-thermal  ionization  must  be  developed 
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and  considered.  An  electron  temperature  of  0.5  volts  would  seem  reasonable  and  results 
in  a  conductivity  of  10  mho/m.  Based  on  this  conductivity  Table  I  below  shows  the  other 
conditions  and  to  calculate  duct  length  as  a  function  of  varying  magnetic  fields  and  for 
three  different  inlet  velocities  corresponding  to  Mach  numbers  of  0.9,  1.275,  2.0. 

Tables  II,  III,  and  IV,  show  the  conditions  and  parameters  which  are  unique  to  each  case. 


TABLE  I 


Mol 

=  17.8 

T,  =  2000°  K 

Pj  =  9. 2  atm 

C 

P 

=  3,893 

joules 

kg°K 

R 

=  466.8 

joules 

kg°K 

=  0.99614 

m° 

r 

=  1.14 

K  =  0.75 

t 

=  0.8 

o 

II 

C 

% 

=  0.2 

TABLE  II 
M  =  2.0 


ui 

= 

2,063  m/sec 

■n'  - 

0.70093 

T 

= 

2,547°  K 

fi  = 

11.617 

°1 

P 

_ 

—  rm  Ift9  watts,  2 

20. 377  x  10  /m 

w  = 

5.3255 

W1 

A1 

= 

2. 4537  x  10"4  m2 

V'  = 

9.860 

Rey1 

= 

1.0452  x  106 

i  = 

0.54011 

fl 

= 

2.90  x  10"3 

q  = 

0. 21053 

r, 

= 

0.1453 
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TABLE  III 


M  =  1.275 


U1 

= 

1,3153  m/sec 

V 

=  0.70093 

T 

_ 

2,  222°  K 

A 

=  11.617 

°1 

P 

W1 

A1 

= 

ii  ha'  ia^  watts,  2 
11.334  x  10  /m 

r 

W 

=  4.649 

= 

4.4116  x  10"4  m2 

V 

=  9.0034 

fl 

= 

1.0452 x  106 

1 

=  0.5163 

= 

2.90 x  10‘3 

q 

=  0.21053 

*1 

= 

-  0. 14533 

TABLE  IV 
M  =  0.9. 


°1 

= 

928.4  m/sec 

n 

T. 

— 

21 1 1°  K 

e 

°1 

P 

7.599 x  109 

w 

W1 

A1 

= 

6.579  x  10-4 

V 

Rey, 

= 

0.7698 x  106 

1 

fi 

= 

3.06 x  10"3 

q  . 

h 

= 

0.0200 

=  0.7395 

=  11.011 
=  4.416 

=  8.7089 

.=  0.50712 


0.21053 
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2  2 

These  were  calculated  for  magnetic  fields  ranging  from  4  webers/m  to  32  webers/m 

which  covers  the  full  range  of  present  or  future  operating  areas.  The  higher  magnetic  field 

is  definitely  out  of  reach  of  any  known  superconducting  material,  however,  it  was  used  as  a 

limiting  point  in  the  analysis.  Figures  9  and  10  show  the  length  and  isentropic  efficiencies 

as  a  function  of  magnetic  field  for  the  three  cases.  The  advantages  of  high  magnetic  fields 

are  obvious,  however  the  increase  both  in  length  and  efficiency  becomes  smaller  per  unit 

2 

increase  in  magnetic  field  above  14  webers/m  .  Therefore  in  the  optimization  based  upon 
weight  considerations  it  appears  the  optimum  value  of  magnetic  field  will  be  near  this  value. 

A  typical  set  of  inlet  and  outlet  parameters  have  been  calculated  and  are  shown  as 
Table  V  to  illustrate  currents,  voltages,  and  thermodynamic  properties. 


M  =  1.275 


TABLE  V 


x,  2. 86  meters 


B  =  12.5 


webers 

2 

m 


T1 

= 

2000°  K 

T2 

= 

1600 

P1 

ty 

= 

9. 2  atm 

P? 

= 

0.69  atm 

= 

0. 99614  kg/m3 

P  2 

= 

9.  295  x  10-2  kg/m3 

ui 

Ai 

i 

= 

1,315  m/sec 

= 

1,  176  m/sec 

= 

4.41  x  10"4  m2 

L 

A? 

= 

49.57  x  10“4  m2 

= 

4. 1 1  x  104  amp/m2 

Z 

i 

= 

4  2 

1.01  x  10  amp/m 

yi 

E 

_ 

1. 23  x  104  volt/m 

y2 

E 

_ 

1 . 10  x  104  volt/m 

yl 

E 

_ 

4.52  x  102  volt/m 

y2 

E 

_ 

50.53  x  102  volt/m 

X1 

x2 

'n 

= 

0.110 

@2 

= 

1.375 

= 

10  mho/m 

n 

= 

2.74  mho/m 
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E.  Effects  of  Strong  Magnetic  Field  on  the  MHD  Process 

A  strong  magnetic  field  is  expected  to  affect  the  MHD  process  in  four  respects: 

1.  The  , "Hal I"  effect. 

2.  Tensor  conductivity. 

®  3.  Hydrodynamic  effects  (wall  friction  and  heat  transfer  coefficients; 

turbulence,  stability,  and  transition  from  laminor  to  turbulent  flow). 

4.  Ion  si  ip. 

These  items  were  discussed  in  our  original  proposal  to  PR -141 108  dated  October 
25,  1963,  which  led  to  the  present  contract.  In  that  discussion  we  cited  31  references 
and  indicated  that  "Hall"  effect  and  tensor  conductivity  were  well  enough  understood 
and,  therefore,  no  further  work  was  proposed.  Hydrodynamic  effects  have  been  studied 
extensively  as  evidenced  by  the  large  number  of  publications  available  in  the  literature. 

In  spite  of  massive  effort,  our  knowledge  is  still  not  as  complete  as  one  could  wish  for. 

In  particular,  information  is  relatively  scarce  regarding  magnetic  effects  on  stability, 
turbulent  boundary  layer,  and  transition  from  laminar  to  turbulent  flow.  This  is  primarily 
a  reflection  of  the  state  of  fluid  dynamics.  Even  without  magnetic  effects,  problems  in  the 
areas  .mentioned  have  to  be  solved  more  or  less  empirically  to  meet  the  needs  of  individual 
situations.  In  our  judgement,  large  research  efforts  spent  on  hydrodynamic  effects  might 
yield  only  small  returns,  if  at  all.  For  this  reason,  we  did  not  propose  to  spend  much  time 
working  in  this  area. 

The  remaining  area,  "Ion  Slip",  is  where  we  shall  concentrate  our  main  effort.  Very 
little  published  information  is  available  on  this  subject.  Nevertheless,  the  effects  of  ion 
slip'  are  believed  to  be  important  when  magnetic  fields  become  very  strong.  This  problem 
should  be  understood  if  strong  fields  are  to  be  effectively  utilized  in  MHD  machines. 

Since  electron  mobility  in  an  MHD  gas  is  much  larger  than  ion  mobility,  the  electric 
current  produced  in  an  MHD  generator  under  a  weak  magnetic  field  is  almost  entirely 
electronic.  As  the  field  increases  in  strength  there  is  a  gradual  transition  from  electron 
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current  to  ion  current.  In  the  limit,  as  the  field  becomes  infinitely  strong  there  can  be  no 
current  whatever.  All  electrons  and  ions  spin  in  infinitesimal  circles.  They  ore  effectively 
arrested.  As  the  magnetic  field  increases  from  zero  to  infinity,  we  expect  that  the  power 
output  of  an  MHD  generator  will  first  increase,  then  level  off  at  a  plateau,  then  decrease, 
and  finally  become  zero. 

in  the  usual  arrangement  of  MHD  generators,  an  ionized  gas  flows  from  a  non-magneti c 
region  to  a  magnetic  region.  As  the  gas  crosses  the  magnetic  field,  the  continuity  of  ion 
flow  requires  that  the  ion  concentration  must  change.  This  situation  is  easy  to  understand  if 
one  considers  the  limiting  case  of  an  infinitely  strong  field.  As  soon  as  ions  enter  the 
magnetic  region,  they  are  forced  to  spin  in  infinitesimal  circles.  Thus,  while  the  neutrals 
are  free  to  move  through  the  magnetic  region,  the  ions  are  arrested.  Therfore,  the  ion  con¬ 
centration  inside  the  field  cannot  be  the  same  as  that  outside. 


Continuity  of  Ion  Flow 


9 
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Consider  the  flow  of  an  ionized  gas  through  a  duct  with  an  electric  field  as 
shown.  Let  the  ion  density  be  n^  and  the  gas  velocity  beu  outside  the  field.  The  ion 
velocity  is  assumed  equal  to  the  gas  velocity  due  to  frequent  collisions.  If  the  ion 
velocity  inside  the  field  region  is  u.  ,  then  continuity  of  ion  flow  requires 

n  u  =  n  u.  (55) 

o  '  i 

where  n  is  ion  density  in  the  field  region.  If  the  ion  mobility  is  p.  ,  then  by 
definition, 


u.  =  u  +  u.  E  (56) 

i  ri  x 

combining  (55)  and  (56)  gives: 


n 

o 


n 


u 


(57) 


Ordinarily  the  quantity  p.  is  small  compared  to  u;  and  the  difference  between  n  and 
n  is  insignificant. 


© 


B 


E 

x 


n 

o 

n 

Extent  of 

r-  u.  ~ * 

Fields 


Continuity  of  Ion  Flow  with  Magnetic  Field 
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If  in  addition  to  the  electric  field  there  is  also  a  magnetic  field  B  ,  then,  as 
will  be  shown  later,  the  ion  drift  velocity  becomes: 


“*t 


I 


(58) 


where  0  .  g  p.  B.  Since  continuity  of  ion  flow  requires 

n  u  =  n  p.  =  n  (u  +  v.  )  ,  (59) 

o  ri  ix  '  ' 

we  obtain  from  (58)  and  (59) 


n 

o 


n 


1  + 


1+*2 


(60) 


the  magnetic  field  becomes  very  strong,  B  becomes  very  large.  E  /uB-e  0  and  . 

2  x 

»  1.  Under  those  conditions  /n  —►0.  Hence  in  a  very  strong  magnetic  field 

the  ion  density  could  be  much  larger  than  that  outside  the  field  region. 


1.  Hall  Effect  on  Ions 

Consider  a  moving  gas  containing  singly  charged  ions.  Let  E  1  be  the  electric 
field  with  respect  to  the  moving  gas;  B  ,  the  magnetic  flux  density;  and  "v".  ,  the  ion 
drift  velocity.  Then  the  force  acting  on  each  ion  (charge  +  e)  is 

F.  =  e  (  I7  +  ~  x  ?)  (61) 

1  1 

If  the  ion  mobility  is  p,  ,  the  drift  velocity  v.  is  given  by 

Mi  F.  =  p.  (  E*  +  v.  xB) 

-  1  ri  1  ' 


v.  = 

1 


e 


(62) 
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j.  =  n.  e  v. 

'■  i  if 


x 


B  acts  in  the  z  -  direction,  then 

(63) 

(64) 

(65) 

(66) 


(67) 
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we  have 


ni  e  Mi 


'IX 


- , —  <  E*  +  0-  E'  ) 

i  +/3.2  x  1  y 


(68) 


n.  e  u. 

j.  = — ' - yi-  (E*  -  0. 

•y  iTT/  y  1 


E;> 


(69) 


By  following  through  a  similar  procedure,  it  may  be  shown  that  the  electron 
current  components  are: 


'ex 


n  e  p 

_i -  •  ( E*  -  P  E'  ) 

i  +  p  2  x  e  y 


(70) 


n  e  li 
e  re 


*y  i  +  4 


0-  ( E*  +  8  E  1  ) 

2  y  ^e  x 


(71) 


where  p  is  the  electron  mobility  and  8  =  p  B. 

Combining  equations  (68)  with  (70)^and  (69)  with  (71),  we  obtain  the  total 
electric  currents: 


x  ix  ex  +  i+  71 


P\  O'*  Pi 


i77  i  +  p2  /  y 


E'  (72) 


i  -  I.  +  i  +  -2i_J  E*  E '  (73) 

y  v  «y  \TT7/  i+  e2}  r  w?  TTT71  x 

i  e/  \  i  e 
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where  (7t  =  n.  e  p.  is  the  ionic  conductivity  and  0 ^ 
electronic  conductivity.  Let 

<T\  (Te 


n  e  u, 

e  re 


is  the 


M  s 


1  +0.2  •  1  +  P ' 


N  B 


CT\  ?\ 

7777 


C7“e  ^e 


1  + 


then  we  have 


i  =  ME'  '+  N  E  ' 

x  x  y 


i  =  ME1  -  N  E ' 

y  y  x 


(74) 

(75) 


(76) 

(77) 


These  equations  are  general.  In  the  work  to  follow  we  shall  apply  them  to  the 
special  cases  of  "Hall"  type  and  segmented -el  ectrode  type  MHD  generators. 


F.  The  Effects  of  a  Non-Uniform  Magnetic  Field 

The  magnetic  field  created  by  a  coil  winding,  whether  made  of  normal  conductors 
or  superconductors,  cannot  be  absolutely  uniform.  In  general,  the  cost  of  the  magnet 
increases  sharply  as  one  strives  for  a  high  degree  of  uniformity.  On  the  other  hand,  a 
nonl-uniform  field  penalizes  the  performance  of  an  MHD  generator.  The  practical 
problem  facing  a  designer,  therefore,  is  to  balance  the  magnet  cost  on  the  one  hand  and 
the  generator  performance  on  the  other  hand.  In  this  study  we  attempt  to  get  some 
quantitative  information  to  help  the  designer  make  his  judgement. 

1.  Past  Work 

The  problem  of  flux  non-uniformity  is  of  course  an  old  one.  ^hat’ particularly 
concerns  us  here,  however,  is  the  effect  of  non-uniform  flux  on  the  MHD  process. 
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(9 ) 

Sonnerup'  '  has  treated  theoretically  the  steady,  laminar  flow  of  a  viscous, 
compressible,  and  slightly  ionized  gas  through  a  magnetic  field.  Although  he 
allows  an  arbitrary  variation  of  the  field  in  the  flow  direction,  his  formal 
solution  ends  with  an  expression  for  the  gas  temperature.  He  does  not 
specifically  discuss  the  effects  of  a  non-uniform  field.  Ehlers^  ^  analyzed  the 
special  case  of  axisymmetric  flow  of  a  conducting  gas  through  a  cylindrical 
channel,  where  the  non-uniformity  in  magnetic  field  arises  from  a  tangentially 
circulating  current.  Sakurai  and  Naito^  ^  ^  investigated  a  more  general  case 
where  a  non-uniform  field  is  generated  by  a  system  of  electric  currents  and  ® 

magnetic  dipoles.  They  obtained  solutions  for  the  variation  of  velocity  and 
pressure  inside  a  two-dimensional  channel,  but  their  solutions  are  restricted  to 
incompressible  flow. 

•  The  effect  of  a  non-uniform  magnetic  field  is  to  some  extent  similar  to 

(12) 

that  of  a  non-uniform  velocity  profile.  In  this  respect  the  work  of  Yeh  and  Sutton'  7 
is  of  interest. 

Although  in  the  present  study  we  are  primarily  interested  in  magnetic  fields 
which  are  non-uniform  in  space,  a  brief  mention  should  also  be  made  of  fields 
which  vary  with  time.  Some  fundamental  work  on  the  motion  of  a  charged  particle 

(13) 

in  a  slowly  increasing  magnetic  field  has  been  done  by  Brder'  .  The  application 

of  a  traveling  magnetic  field  tc  accelerate  a  plasma  has  been  studied  fairly 

(14)  (15) 

extensively.  Some  representative  publications  are  those  of  Meyer  7,  Light'  , 
and  Smotrich,  Janes,  and  Bratenahl^^. 

2.  Model  and  Approach 

The  question  that  we  are  particularly  interested  in  here  is  the  effect  on  MHD 
generator  performance  of  a  non-uniform  magnetic  field.  The  literature  already 
cited,  although  helpful,  is  not  directly  applicable.  Before  a  meaningful  analysis 
can  be  made,  let  us  first  examine  the  nature  of  non- uniformity  produced  by  a  typical 
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electro-magnet.  Figures  11,  12,  and  13,  taken  from  reference  (17),  show  the  flux 
density  variation  throughout  the  working  volume  of  a  rectangular  MHD  duct. 

Figure  11  shows  the  variation  in  the  flow  direction.  Figures  12  and  13  show, 
respectively,  the  variations  in  the  horizontal  and  vertical  directions  (normal  to 
flow).  An  examination  of  these  figures  reveals  that  the  largest  variation  occurs  in 
the  flow  direction.  Furthermore,  almost  all  of  this  variation  takes  place  in  the  first 
quarter  and  last  quarter  of  the  duct  length.  ^The  field  is  quite  uniform  in  the  center 
portion  of  the  duct.  For  the  purpose  of  our  analysis  to  follow,  therefore,  we  will 
adopt  a  model  as  shown  in  Figure  14.  Between  0.  and  L/4,  the  flux  density 
increases 'from  rfB  to  B  ,  where  a'is  some  number  less  than  unity.  Between  L/4 
and  (3/4)  L,  B  =  Bq  =  constant.  Between  (3/4)  L  and  L,  B  decreases  from  B^  to 
a  Bq  .  From  0  to  L/4  the  variation  of  B  may  be  represented  by  the  equation: 


1  -  a' 
TS3TT 


a7-  .9689 

.0311 


(78) 


Equation  (78)  will  be  used  in  our  analysis  to  be  reported  next  quarter.  We 
plan  to  vary  the  degree  of  non-uniformity  by  changing  the  number  a'.  If  a'=  1, 
equation  (78)  reduces  to  B/B^  =  1. 
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III.  SUPERCONDUCTING  MAGNET 

Refrigeration  Requirements  for  Various  Optimized  Methods  of 
Cooling  Cryogenic  Electrical  Current  Leads,  Conduction  and 
Radiation  Transfer  Paths.  (W.  A.  Stewart) 


A.  INTRODUCTION 

In  cryogenic  systems,  and  in  particular  the  superconducting  magnet  cryogenic 
electrical  systems,  heat  is  transferred  from  an  ambient  high  temperature  region,  through 
vacuum  insulation  by  radiation  and  conduction  in  support  structure  and  is  transferred  by 
conduction  and  generated  by  Joule  heat  loss  in  electrical  leads  carrying  current.  A 
superconductor  must  be  maintained  near  4.2°  K,  the  standard  vaporization  temperature 
for  helium.  Although  the  usual  cryogenic  techniques  result  in  rather  small  heat  fluxes, 
because  refrigerator  coefficients  of  performance  are  small  (  <  Tq  /Tq  -  T^  ,  where  a 
refers  to  rejection  and  o  refrigeration  temperatures)  a  much  larger  refrigeration  work 
rate  or  power,  W  ,  is  required.  Nearly  1000  watt  of  heat  at  4.  2°  K.  Heat  that  is  ^ 
rejected  therefore  is  practically  equal  to  W. 

The  particular  superconducting  magnet  system  under  study  is  to  provide  a  high 
field  across  a  long,  large,  varying  section  MHD  duct  which  is  part  of  a  nuclear-MHD 
space  power  system.  The  one  important  subsystem  comprises  (1)  MHD  duct  insulation, 

(2)  cryogenic  insulation,  (3)  magnet  winding,  (4)  magnet  force  supports,  (5)  refrigeration 
apparatus  and  (6)  heat  rejection  radiators.  The  weights  of  refrigeration  apparatus,  the 
penalty  to  weight  of  the  power  used  for  refrigeration  and  heat  rejection  radiators  are 
proportional  to  W.  Refr  igeration  power  studies  are  reported  here,  for  various  configura¬ 
tions  of  cooling  which  will  minimize  that  power.  Although  the  analytic  work  presented 
here  does  not  relate  to  the  details  of  constructed  systems  which  will  accomplish  the  many 
different  methods,  separate  studies  aie  being  made  to  find  the  relative  practicality  and 
details  for  implementing  the  theoretical  methods  derived  in  this  work.  It  is  sufficient  to 
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$ay  that  the  methods  can  be  used  in  ways  that  are  now  practised  by  we  or  others  in 
cryogenic  systems  and  modifications  are  feasible.  The  purpose  of  the  study,  moreover, 
is  to  obtain  comparisons  of  methods  which  will  lead  to  best  or  optimum  (for  weight  and 
reliability)  magnet  systems. 

It  is  assumed  that  an  MMD  magnet  may  require  provision  for  current  charging  in 
the  space  application.  Therefore  facility  for  cooling  electrical  leads  which  deliver 
the  large  amount  of  magnetic  energy  will  be  required.  This  is  notwithstanding  the  fact 
that  persistent  mode  operation  eliminates  the  need  for  leads  when  steady  state  is  attained 
after  a  long  time  for  charging. 

B.  TECHNICAL  DISCUSSION 

Cooling  of  (1)  electrical  leads,  (2)  conduction  paths  and  (3)  radiation  paths  are 
studied  analytically  and  in  detail  for  varying  or  constant  electrical  resistivity,  thermal 
conductivity  and  ratio  of  actual  COP  (coefficient  of  refrigerator  performance)  to  Carnot 
COP,  ^  ,  in  part  C  •  When  some  heat  flux  is  removed  at  a  temperature  above  the  lowest 
temperature,  Tq  ,  refrigeration  has  a  higher  COP  and  total  work,  W,  is  reduced.  This 
technique  is  exploited  by  (a)  optimum  staged  cooling  all  along  a  heat  transfer  path, 

(b)  vapor  cooling  to  remove  the  same  heat  in  each  temperature  increment  along  a  heat 
transfer  path  and  (c)  two  stage  cooling  at  an  intermediate  temperature,  T^  ,  between  the 
high,  Tq  ,  and  low,  Tq  .  In  addition  (d)  cooling  can  be  applied  just  at  the  cold  end, 

T  .  All  of  the  methods,  a,  b,  c,  d  are  applied  to  path  types  1,  2,  3  in  the  analysis 
and  further  optimization  of  length-area  ratio,  I/A,  Tj  ,  etc.  found. 

1 .  Electrical  Leads 

McFee,  in  reference  18,  examined  cooling  of  an  electrical  lead  at  Tq  and  found 
optimum  I/A.  Part  C,  Section  3,  Case  A  essentially  reproduces  this  result.  Kothmann 
and  Stewart,  in  reference  19,  examined  vapor  cooling  and  the  results  are  contained  in 
Part  C,  Section  2,  Case  A.  Sobol  and  McNichol,  reference  21,  later  gave  experimental 
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verification  of  their  analysis  of  this  cose.  Tsu,  in  reference  22,  examined  optimum 
cooling  but  did  not  analyze  by  a  variational  approach.  He  achieved  about  a  40% 
improvement  over  the  McFee  result. 

The  results  of  the  present  analysis  considering  the  model  electrical  lead  shown 
in  Figure  15,  are  shown  in  Figure  17  for  constant  properties,  and  in  Figure  18  for 
properly  varying  resistivity.  It  may  be  noted  in  Figure  18  that  the  optimum  case  is 
really  slightly  poorer  than  for  vapor  cooling.  This  is  the  result  of  the  limiting 
restriction  in  the  problem  where  calculus  of  variations  is  used  that,  a  given  heat  flux 
or,  afternatively,  T'  (o)  is  taken.  The  parameter  f  relates  heat  flux  at  the  bottom 
of  the  lead  to  the  cooling  along  the  lead.  Figure  17  shows  the  strong  influence  of 
this  parameter. 

Smith,  in  reference  20,  used  McFee1  s  approach  and  found  that  for  varying 
properties  and  cooling  at  T  ,  the  minimum  refrigeration  work  would  be  40  watts  per 
ampere  in  the  lead  if  COP  is  0.001  from  4.  2  to  300°  K.  Figure  18  shows  that  a  10.7 
reduction  is  possible  so  that  only  3.64  watts  per  ampere  might  be  required. 

2.  Conduction  and  Radiation  Paths 

Analysis  of  cooling  of  conduction  and  radiation  paths  is  contained  in  Part  C, 
Sections  4  through  8.  Results  are  shown  ploted  in  Figures  19  and  24.  Again,  as  seen 
in  Figures  19,  21,  and  24  for  small  values  of  T ^/T  ,  the  restriction  imposed  by  the 
natural  boundary  condition  of  the  variational  problem  makes  the  "optimum"  work  more 
than  work  with  vapor  cooling. 

No  complete  conclusions  will  be  drawn  except  those  that  are  self  evident  from 
the  curves.  Substantial  reductions  in  work  and  system  wejght  are  possible  when  using 
an  optimum  method  rather  than  cooling  only  at  T  .  Optimum  two  stage  cooling  is  nearly 
as  effective  as  vapor  cooling.  In  general,  vapor  cooling  is  also  near  optimum  cooling 
except  for  radiation  exchange.  With  radiation,  the  high  order  transfer  coefficient 
variation  with  temperature,  makes  optimum  cooling  much  better  than  any  other  method. 


C.  ANALYTICAL  DISCUSSION 
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1.  Optimum  Cooling  of  on  Electrical  Conductor 

One  dimensional  heat  flow  in  a  conductor  such  as  shown  in  figure  16,  with 
current,  i,  electrical  resistivity,  f ,  and  thermal  conductivity,  k,  is  governed 
by  equation  79  . 


0 


(79) 


The  heat,  d  x,  is  removed  from  the  length  element,  dx,  at  temperature  T  by 
refrigeration  and  is  rejected  at  the  ambient  or  radiation  temperature,  .  Therefore 
refrigeration  work  for  ^  ,  using  equation  79  is: 


The  variational  problem  to  find  T  and  hence  ^  which  minimizes  W  has  the  following 
form,  Euler  equation  and  natural  boundary  conditions: 


F  is  the  integrand  of  equation  80. 
T  is  to  be  fixed  at  x  =  0  and  I  , 


Both  natural  boundary  conditions  are  satisfied  since 
3F/£Tm  is  equal 


-  1)  kA  =  0  at  x  =  I  , 


2* 


T  =  Tq  and  dT/dx  or  heat  flux  may  be  considered  fixed  at  x  =  0.  When  \  is  constant, 
the  Euler  equation  becomes: 


(kA)  T 

■W" 

T2 

Tnm 


E 


i2  P 

c)(kA) 

“Jr 


+  T1 


,2  waJ1=  0 


(82) 


The  partial  differentiations  are  taken  when  kA  or  f / A  are  explicit  functions  of  T.  Most 
materials  of  interest  as  electrical  leads  (having  a  small  P k  product)  also  have  nearly 

constant  thermal  conductivity.  Then  introducing  dimensionless  temperature,  O  =  T/T  , 

2  a 

and  parameter  y  =  i  f^/  2  (kA)  ATq  ,  while  A  cannot  be  an  explicit  function  of  T, 

reduces  equation  32  to: 
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,2 

0 
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-  (0  e2)  r 


MW  fig)  _  o 

do 


(83) 


Two  cases  are  of  interest  and  for  which  the  non-linear  Euler  equation  possesses  analytical 
solutions.  These  are  for  resistivity  constant  and  for  resistivity  varying  linearly  with 
absolute  temperature. 


Case  A  -  Resistivity  Constant 

The  equation  and  boundary  conditions  are: 

e"  =  ®ll  e(0)  =  eo,e(i)  =1 

0 

Reference  23  gives  the  first  integration  of  this  equation  as: 

2  +  2 

©'  =  2y©  -  r  a  O 


(84) 


(85) 
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Either  the  positive  or  negative  sign  may  be  taken,  depending  on  the  constant  of 
integration  which  is  y  a.  When  the  positive  sign  is  taken  final  integration  and  the 
boundary  conditions  give: 


(86) 

(87) 

The  total  refrigeration  work,  including  that  due  to  the  heat  flux  (  k  A)  T'  (0)  and 
is  found  as: 


ycTy'x  =  sin  ^  (  a  0  -  1)  -  sin  ^  (a  0  -  1) 

o 

/a  y'  I  =  sin  ^  (  a  -  1)  -  sin  ^  (  a  -  1) 


W  = 


kAT 


1C* 


D(©  +2y) 


dx 

60 


d  0 


o 


0  (o) 


} 


(a-1)  -  sin 


(88) 


Here  r  is  the  ratio  T  /T  =  1/ Q  .  The  optimum  value  of  I  for  minimum  W  can 
be  found  alternatively  by  setting  the  derivative  of  equation  88  with  respect  to  a  equal 
to  zero  and  solving  for  a.  That  optimum  a  is  2,  and  optimum  |  is  found  from  87. 
Then: 


4  /T77  -  7  f 


2  sin 


(89) 
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When  the  negative  sign  in  equation  85  is  used  the  solutions  corresponding  to  equations  86,  88, 
89  ,  90  ,  and  91  are: 


x  =  In 


^  a  ©  +  a^  0  ^  +  a  ®  +  1 


/  2  Z 

2a9  +  a  0  +a©  +1 

o  o  o 


(86  a) 


(88a) 
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Case  B  -  Varying  Resistivity 


f~  />»(T/T0) 


Equation  83  becomes  now: 

,2 


0"  =  -rQ  ~  ~  Y  0  >  ©(»)  -  /  ©(  I  )  -  1 


Reference  23  gives  the  first  integration  of  this  equation  as: 

2 


*'2  -  ^  * 


2  y  ©‘ 


(92) 


(93) 


The  constant  of  integration  is  2  y  //S  and  is  positive.  Integration  of  equation  93 
and  the  boundary  conditions  give: 


W1- ft  •l) 

(i  +y  i  -  to' ) 


( i  Vi 

(  1  +^1  ‘  )_ 


(94) 


(95) 


The  total  refrigeration  work,  found  in  the  same  manner  as  in  equation  88  is: 


Optimizing  W  to  find  /3  ,  which  again  is  equivalent  to  finding  optimum  I,  results  in  an 
optimum  (3  of  1. 


W 


opt. 


(97) 
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2.  Vapor  Cooling  of  an  Electrical  Conductor 

When  a  lead  is  cooled  by  a  stream  of  vapor  flowing  along  the  conductor  from  cold 

to  warm  end  and  convection  heat  transfer  coefficients  are  high,  then  vapor  temperature 

and  lead  temperatures  are  the  same  at  any  location,  x.  Furthermore,  since  vapor  specific 

,  * 

heat  C  is  independent  of  temperature,  (p  ,  is  constant.  That  means  that  the  same 
cooling  or  heat  removal  occurs  between  equal  temperature  increments  all  along  the  lead. 
Since  d>  (heat  removed  per  unit  length)  is  equal  0  (heat  removed  per  unit  temperature 
change)  times  d  T/d  x,  equation  79  may  be  written  as: 

TC"  (kA'rf')  +-JT  "  t  -ITT  =  0  (100) 

Again  two  solutions  are  of  interest,  either  constant  or  linearly  varying  resistivity  and 
constant  conductivity,  k,  and  area  A. 

Case  A  -  Resistivity  Constant 

i* 

When  9  =  kAD,  a  constant,  equation  100  and  boundary  conditions  become: 

y  /  ©(o)  =  ©Q  ,  ©  ( I  )  = 


=  D&'  =  -  2 


1 


001) 
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The  solution  of  this  equation  is: 


0-eo-^;)adx  -(!-©;* 


©  = 


•  Dl 
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“di  r 

e  -  1 


TT 


002) 


Total  refrigeration  work,  including  that  due  to  a  heat  flux  at  the  cold  end  of  the  lead  is: 


(103) 


Using  equation  103,  W  may  be  minimized  with  respect  to 
for  optimum  I  is: 


+ 


1  -  e 


-Dl 


=  0 


The  resulting  equation 


(104) 


Assuming  that  D  I 


is  large,  e  is  negligible  compared  with  one.  Then 


Dl 


opt. 


(105) 


When  vapor  cooling  is  used,  there  is  often  a  unique  relation  between  the  heat  transferred 
to  evaporating  refr  igerant  and  the  vapor  cooling  parameter  D.  For  a  flow  rate,  m. 
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caused  by  heat  flux  kAT1  (o): 


m  Cp 

Ea 


loCp  ©‘(a)  = 

hfg  ©o 


(kAT' (o)  /  hfg)  Cp 
- &K  - 

©•  (a) 


(106) 


Here  f  is  a  factor,  near  1  for  helium,  defined  as  heat  of  vaporization,  h^g  ,  divided 
by  C  Tq  .  An  equation  using  the  optimum  result,  equation  105^  and  equation  102 
in  equation  106gi''es: 


(107) 


2 

Equation  107  is  satisfied  for  D  /  2  y  nearly  equal  r/f  ,  which  is  generally  larger  than 
one.  Then,  equation  105  shows  that  optimum  D  I  is  larger  than  two.  This  conclusion 
justifies  neglecting  e  ^  in  equation  104, 


The  results  may  finally  be  written  making  use  of  equations  106  and  107  in  equations  103 
and  105. 


(108) 


(109) 

(110) 
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where  g  =  1 


-e  -  (  1  +  JJ—)  * 


An  optimum  f,  equivalent  to  providing  refrigeration  at  just  Tq  and  also  vapor  cooling 
so  that  f  is  variable,  can  be  found  from  equation  108.  When  g  equals  1. 


Case  B  -  Resistivity  Variable,  ft=  f*  (T/Tq) 

When  (f)  =  2kAd,  a  constant,  equation  100  and  boundary  conditions  become: 

e"  -  2dQ'  +  2  y©  =  0,  ©  (o)  =  ©q,  ©  ( I  )  =  1  I 

The  solution  of  equation  1 13  is: 


©dxf_  ax  .a  -d  .  smh  ax 
=  e  ©  e  -(©e  -e  )  — s-r — r- 
L  O  '^o  '  Sinhal  _J 

=  e^^©^  sin  ft  (  I  -  x)  +  e  ^  sin  xj 


sin  p  I 


where:  a  =  d  -  2  y  ,  ft  =-J  2  y  -  d* 


The  total  refrigeration  work,  considering  that  d  =  ©  (o)  /  2  f  0  ,  is: 

W  =  2  k  A  To  |*|n  r- (1  .  f)  (  1  .  _L)|  d 
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To  minimize  W  it  is  necessary  to  minimize  d  with  respect  to  I.  Using  the  equation  for  d  as: 

a 

1  e 

o 

sinh  a  I 


,.©'(»)  deo+a 

°  7T5 - 
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,  -  al  -d  I  ' 
( ©  e  -  e  ) 
o 


(116) 


From  equation  116  these  alternative:  conditions  are  found: 

■(dl)  r  1 


(2 .  -  1)  -  _(1  -if 
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r.e  v  -  cosh 
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sinh  1 - (d  I)  | 
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(117) 


(2  f -  1)  = 


-  1 


-(dl) 

re  -  cos 


iW7 


1  (dl) 


sin (  — -  1  (d  I) 


(118) 


j2 . 


If  d  is  greater  than  but  near  2  y,  equation  117  reduces  to: 

=  T3TT-  [re'<dl)  -  ’] 


(119) 


Since  it  is  possible  to  satisfy  equation  119,  and  solve  it  for  (dl),  d  may  be  as  small  as 
2y.  The  second  condition,  equation  1 18,  therefore  should  be  used  to  find  dl  which 
minimizes  d  or  maximizes  2  y/d2  in  equation  118.  An  alternative  is  to  find  maximum 

I  ^  1 

-J(  2  y/d  )  -  1  =  p  optimizing  d  I.  Equation  118  can  be  solved  for  r  and  the  variation 
of  it  set  equal  to  zero  for  p  and  d  I  varying.  Then  dp/d  (d  I)  is  set  equal  zero  from  that 
expression  to  find  an  optimizing  condition.  This  results  in  an  equation  for  pdl  in  terms 
of  p.  When  f  =  1,  dl  =  I  nr  and  p  may  be  found  from  the  original  equation  118  by  a 
graphical  solution  of  the  transcendental  function. 
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3.  Cooling  an  Electrical  Lead  at~One  End 

An  electrical  lead  that  is  cooled  only  at  x  equal  zero  is  described  by  equation 79 
except  that  the  advantages  in  staged  cooling  are  not  used  and  if)  is  zero.  Again  two 
cases  are  of  interest,  constant  or  linearly  varying  resistivity  with  constant  conductivity 
and  area. 


Case  A  -  Resistivity  Constant 

Equation  79  and  its  boundary  conditior^  become: 

©"  +2y  =  0  ,  ©(c)  =  0o  ,  ©(I)  =  1 


The  solution  of  equation  120  is: 


©  = 


©  + 
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a  -  ©0) 


Y  1 


(120) 


(121) 


The  refrigeration  work  required  to  remove  the  flux  of  heat  at  x  =  0,  Tq  lo  Tq  is: 


,  .  ,  T  -  T 

W-  ±±(  \  °l  T'(o> 


n 


(122) 


To  minimize  W  it  is  only  necessary  to  minimize  ©'  (o)  with  respect  to  I. 


©'  (o)  = 


1  “ 


+  yl 


(123) 


From  this  the  optimum  length  to  area  ratio  and  refrigeration  work  are: 


\ 


11  „  I  2  k  (T  -  T  ) 1  /  2k  T  '  .  . 

opt  I  'a  0 I  a  I  r  -  1 


(124) 
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Case  B  -  Varying  Resistivity  P  =  PQ  {T/Ta  } 

Equation  79  and  its  boundary  conditions  become: 

©"  +  2y©  =  0,  ©(o)  =  0Q  ,  0(1)=  1 

The  solution  of  equation  126  is: 

0q  sin  ^2  y  (  I  -  x  )  +  sin  «^2y  x 

©  =  ...  — —  i  .  — — — . . - 

sin  ^2  y'  I 


(125) 


(126) 


(127) 


The  refrigeration  work  required,  that  given  by  equation  122  ,  can  be  minimized  with 
respect  to  I.  Then: 


4.  Optimum  Cooling  of  a  Conduction  Path 

One  dimensional  heat  flow  in  a  conduction  path  is  shown  schematically  in  Figure 
The  heat,  fi  d  T,  is  removed  at  temperature  T  by  refrigeration  and  rejected  at  ambient 
Tq  .  Heat  flux  Q  flows  toward  the  cold  end  at  the  section  x  with  temperature  T.  Total 
refrigeration  work  is: 
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Heat  transfer  can  be  described  in  two  ways: 

,  kAdT 

dx  "  “57TT 

(131) 

d  Q=  f  dT 

(132) 

Equation  131  establishes  a  constraint  for  the  minimization  of  W.  The  integral  of  equation  131 
is: 


o 


k  A  d  T 

imrr 


(133) 


Using  equation  132  inequation  130  and  introducing  the  Lagrange  multiplier,  X  ,  the 
problem  is  to  minimize: 


The  integral  has  the  form 


I  = 


f2(T)  TT-  +  x 


f,  0) 


+  c. 


dT 


(135) 


The  Euler  equation  is: 


df2  (T) 


Xf,  (T) 
Q2 


(136) 


dT 


0 
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Consequently : 


q2  =  (137) 

a 

The  Lagrange  multiplier  is  evaluated  by  using  equation  137  in  equation  133  .  The 
result  from  Q  is: 


o 


Refrigeration  work  can  be  found  as: 


(138) 


Two  cases  are  of  interest,  that  for  variable  thermal  conductivity,  k  ,  (actually,  k  A) 
and  variable  efficiency,  ,  and  that  for  constant  k  and  ^  . 

Case  A  -  Variable  Conductivity,  kA  =  k^  A(T/TQ)n,  and  Efficiency,  ^ 

(T/Tq  )m 

When  the  conductivity  and  efficiency  functions  are  used  in  equations  138  and  139 
the  results  are: 
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Solution  of  equation  131  for  for  resulting  temperature  distribution  with  Q  given  by 
equation  138  and  the  property  variations  is: 


Case  B  -  Constant  Conductivity  and  Refrigerator  Efficiency 
Application  of  equations  138  and  139  gives: 


a 


The  resulting  temperature  distribution  is: 

*  is  constant,  vapor  cooling  is  also  optimum  cooling. 


(142) 


(143) 


(144) 


(145) 


Note  that  since 
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5.  Vapor  Cooling  of  a  Conduction  Path 

When  vapor  cooling  is  used  (f)  is  a  constant.  Recalling  equation  106  where  f  was 
introduced,  equation  130  may  be  written 


W  =  <p* 


t(t  _1)  dT  +wr  W 


L  To 


(144) 


Equation  132  may  be  used  to  find  Q  as: 


Q  =  0o  *  f  <T-T0)“  f  [T-(l-f)  Tj 


(147) 


Then  equation  131  when  integrated  evaluates  ^  as: 


f- 


kAdT 


i  [t-o-otJ 


(148) 


When  properties  are  variable,  it  may  be  seen  that  k  variations  appear  only  in  ^  and 
*L  variations  only  in  the  integral  of  equation  146  for  W.  Several  useful  results  may  be 
obtained. 


Case  A  -  Refrigeration  Efficiency  Variable  ^  =  (  T  /  T  ) 


m 


Refrigeration  work  from  equation  146  becomes: 


W  = 


m 


(  m  -  n  + 


m 

r 


im  -  1 


+  f  r 


m  -  1 


(149) 


There  are  three  useful  cases  when  considering  conductivity  variation. 


44 


Case  A-1  -  Conductivity  Variable,  k  =  kQ  (T  /  Tq  )  ,  f  =  1 


This  case  corresponds  to  the  use  of  helium.  Equations  148  and  149  give: 


k  A 

f  =-T- 


—  (  1  -  — - ) 

n  '  n  ' 

—  r  — 


_  a  a  /  l  ,  ,  1  ^  1  ^  m  +  1  m  ^  m 

W=  -= — T —  ( —  (1 - )  — > — tv-  +  - r  + - 

I  in  n  (m-1)  m  m  - 

Case  A-2  -  Conductivity  Variable,  k  =  k^  (T /  Tq  ) 

Integration  in  equation  I48may  be  performed  with  f  general. 


,*  k  A 

t  *-r- 


0-7-)  +  (i 


k  AT 
a  a 


+  (1-  f) 


-(i-o 


1  ,  /  1  ,  r  »  m  ,  1  .  , .  m  -1 

sttsttv  +  <—  *  f> '  -  +  ,)r 


Case  A- 3  -  Conductivity  Constant 

For  this  case  equations  148  and  149  give: 
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p_  kc,ATa  J|  r-  (1-0 ”1 F  1 
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Case  B  -  Refrigeration  Efficiency  Constant 

Refrigeration  work  from  equation  146  becomes: 


Inr-  {1  -  f)  (1 


056) 


Again  there  are  three  subcases  of  interest.  Since  the  same  ^  1  s  as  found  in  subcases 
A-l,  A-2  and  A-3  apply,  only  refrigeration  work  will  be  summarized. 


CaseB-1  -Conductivity  Variable,  k  =  k  (T/T  )  ,  f=  1 

— — —  -  a  a 


w  = 


k  AT 
a  a 


_1_  (  i  -  J_)  In  r 

n  n 

r 


(157) 


Case  B-3  -  Conductivity  Constant 


W  = 


k  AT 


a  a 


(159) 


6.  Cooling  a  Conduction  Path  in  Two  Stages 

A  conduction  path  cooled  in  two  places,  a  the  cold  end  and  somewhere  (x  =  Iq)  at 
temperature  T^  is  an  approximation  of  optimum  cooling.  For  this  problem  both  the  best 
temperature,  T^  ,  and  location  of  intermediate  cooling,  I  ,  may  be  found.  The  general 
case  where  k  and  are  variable  as  well  as  that  special  case  n  =  1  and  ^  constant  are 
of  interest. 
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Gose  A  -  Variable  Conductivity  ,  k  =  kQ  (T/Tq  )  n  and  Efficiency/  %  -  (T/Tm)  111  • 


The  equation  governing  heat  transfer  in  length  I  is  equation  131. 


k  A 

dx  =  — - -  Tn  dT,  T  (o)  =  T  ,  T(I)  =  T 

_  _  n  o 

Q  T 
o  a 


a  a 


(160) 


Integration  gives: 


Q  = 
o 


k  AT  . 

a  a  1 


(n  +  l) 


»  n  +  1  /  j  \  n  +  1 
1  \  /  c 


(161) 


Similarly  if  the  heat  flux  Q,  is  flowing  from  T^  to  Tj  : 


Q1  = 


k  AT 
a  a 


ttht 


i  - 


n  +  1 


(162) 


The  work  required  by  refrigeration  of  O  and  Q.  -  Q  ,  i.e.,  W  and  W,  is: 

O  I  o  0  1 


<V  Wl*4- 

u  a 


TJ--1J(Q]-Q0)  063) 


Total  work  required  is  Wq  +  and  when  r  =  Tq/  Tq  ,  S  =  Tq  /  Tj  and  X  =  I  ^  /  Iq  : 


(164) 
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Then  also: 


(165) 


066) 


The  optimum  value  of  X  is  found  as: 
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2 

opt. 


(s-os™  a-^Vrr) 


(r  .  ,)  r  m  .  (S  -  1)  s  liT-Ur 
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Then  minimum  work  is: 
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k  AT 
a  a 


°pt:  7  a  '  \n+1 


(S-l)Sm  (1  -  1 


n+1 


For  any  n  and  m,  numerical  trial  is  the  most  effective  method  to  find  best  values  of  S 
given  r. 


Case  B  -  Linearly  Varying  Conductivity,  k  =  kQ  (T/Tq),  Efficiency  Constant 


From  equations  167  and  168  for  n  =  1  and  m  =  0. : 


(169) 


070) 
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Case  C  -  Constant  Conductivity  and  Efficiency 

From  equations  167  and  168  for  n  =  0  and  m  =  0 


X 


opt. 


S-  1 
r  -  S 


W 


k  AT 
a 


a 


S  -  1 

■fT 


+ 


r  -  $ 
i frS ' 


2 


071) 

(172) 


The  optimum  value  of  S  is  found  from  equation  172  as  S  =V^  •  The  refrigeration  work 
is: 


W 


k  AT 
a  a 


opt. 


r1/4-  1 


n  2 


17T 


073) 


7.  Cooling  a  Conduction  Path  at  One  End 

The  refrigeration  work  required  for  a  conduction  path  cooled  just  at  Tq  is  a 
specialized  case  in  the  preceding  section  where  S  =  r.  Use  of  equation  168  then  gives: 


(174) 


This  result  may  be  used  where  refrigeration  efficiency  and  thermal  conductivity  vary 
(n  f  0,  m  f  0),  or  where  one  or  both  these  parameters  are  constant  so  that  m  =  C  or 
n  =  0. 


8.  Cooling  a  Multiple  Shielded  Radiation  Path 

Heat  transfer  by  radiation  between  a  surface  No.  1,  and  a  second  surface,  No.  2, 
very  near  to  it  is  given  by: 


Q  = 


V'  A 


<t/-t24>  =(T7n5  a(Vt2> 


175 


Equation  175  defines  an  apparent  conductivity  k^  .  If  there  are  a  number,  N,  of 
radiation  shields  between  Tq  and  Tq  ,  then  adjacent  shields  or  radiation  surfaces  have 
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absolute  temperatures  that  are  not  very  different  in  magnitude.  Because  of  this,  kf  may 
be  written  as: 


»■!  (y  Tjd,2  +  T,2)  =  4_CM  y3 

N  <rnr  - 


076) 


Emissivities,  €  ,  are  made  much  less  than  1,  so  that  a  good  approximation  for  k  is 

3  r 

4  £(T)  (t“  I  T  /N.  Now  it  is  possible  to  use  kr  in  the  preceding  theory  for  conduction 

paths  with  variable  conductivity.  Now  if  €/ £q  =  (T/TQ)  j  ,  n  =  3  +  j.  For  example, 

with  constant  ratio,  n,  of  refrigeration  efficiency  to  carnot  cycle  efficiency  and  constant 

emissivity  €.  . 


Cose  A  -  Optimum  Cooling 

From  equations  140  and  141  with  m  =  0,  n  =  3: 


Case  B  -  Vapor  Cool  i  ng 

From  equations  148  and  156: 


W  = 


-0-0(1- 


(177) 


(178) 


079) 


'(180) 
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Cose  C  -  Two  Stage  Cooling 

From  equation  168,  with  X  taken  as  optimum  and  interpreted  as  /  Nq  ,  the 
number  of  shields  between  hot  Tfl  amd  cooling  at  Tj  divided  by  the  number  of  shields 
between  T,  and  cooling  at  T  : 

1  o  ;  „ 


At  (TAJ 


opt.  -0 


AT  H  IT  T 

( T"))  ( s_1)  (  1  "T^  +(r_s)(-^ — 3-) 

L  S  -l  L  S  r  - 


It  may  be  shown  that  optimum  S  is  found  from: 


4  C4  c4  1 

r  -  5  S  -  I 

~r ;  zr  s  - 1 

r  (r-  S) 

4  4 

When  r  is  larger  than  about  3,  S  «  r  and  equation  182  may  be  solved  for  r  when 


S  is  optimum. 


*  2  4 

S  \  (S*  -  1) 

r  =  S  +  °£  :  ■.  °Pn 

opt.  16  S  -  1) 

r  opt. 


The  resulting  expression  for  W  in  terms  of  S  j,  becomes: 

4*<?-AT  4  f  .  (S4  -  1)  (5$  -  4)  2 

VA/  _ _ a  )  1  opt.  opt. 


S  f  (S  .  -  1) 

opt.  opt. 


opt.  N  \  48 


Case  D  -  Cooling  only  at  Tq 
Equation  174  gives: 


(185) 
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Table  of  Symbols  -  Equations  1  through  78 

Cross-sectional  area  of  duct 
Magnetic  flux  density 

Magnetic  flux  density  in  central  region  of  duct 

Constant  pressure  specific  heat 

-19 

Electronic  charge  (1.6x10  coul) 

Applied  electric  field 

Electric  field  relative  to  gas 

Frictional  coefficient 

Force  acting  on  an  ion,  equation  61 

Convective  heat  transfer  rate 

Current  density 

Boltzman' s  constant  (1.38  x  10  joule/°K) 

Loading  parameter,  equation 9 

Constant,  equation  41 

Constant,  equation  29 

Mach  number 

Constant  in  equations  74  -  77  only 
Number  density  of  particles 
Constant,  equation  28  -  30  only 
Drag  at  wall,  equation  5 
Constant  in  equations  74  -  77  only 
Static  pressure 

Convective  power  (inlet)  per  unit  area,  equation  20 

Constant,  equation  40 
Reynolds  number,  based  on  duct  width 
Hydraulic  diameter 
Duct  perimeter 


t 

T 

T 

°1 

u 

V 

w 

X 

o 

X 

y 


Temperature  ratio,  equation  13 
Static  temperature 

Total  or  stagnation  temperature  (inlet) 

Velocity 
Ion  drift  velocity 
Constant,  equation  25 
Constant,  equation  19 
Distance  along  axis 
Variable  f(  t  ),  equation  38 


a  ' 

e 

0 

e 

r 

e 

n 

n' 

% 

K 

^e 

t 

P 

<7“ 


Frictional  parameter,  equation  11 
Magnetic  field  parameter 
Constant,  equation  14 
Product  (  p^B) 

Product  (p.B) 

Ratio  of  specific  heats 

Seeding  ratio  (moles  cesium/moles  lithium) 

Ratio  of  electrical  power  to  total  inlet  enthalpy,  equation  49 
Constant,  equation  8 

Isentropic  or  turbine  efficiency,  equation  44 

Viscosity 

Electron  mobility 

Ion  mobility 

Constant,  equation  34 

De.isity 

Scalar  electronic  conductivity 
Electronic  conductivity 


Ionic  conductivity 
Constant,  equation  27 
Constant,  after  equation  39 
V  Constant,  equation  35 

Superscripts 

—  Designates  vector  quantities 

Subscripts 
e 

1 

x 

y 

z 
1 

2 

Exponents 

y  Conductivity  variation  with  temperature,  equation  22 

z  Conductivity  variation  with  pressure,  equation  23 


Electron 

Ion 

Coordinate,  axial  direction  of  flow 
Coordinate,  direction  of  applied  electric  field 
Coordinate,  direction  of  applied  magnetic  field 
Specifies  inlet  conditions 
Specified  outlet  conditions 
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Table  of  Symbol>  -  Equation!  79  through  185 

Cross  section  area  of  lead,  conduction  path  or  radiation  path 
Specific  heat  of  vapor 

Cooling  parameter  for  lead  with  constant  resistivity,  =  Q*  /  kA. 

Cooling  parameter  for  lead  with  variable  resistivity,  =  fi*/2  kA. 

Integrand  of  equation  2. 

Parameter,  =  Q  /  T  m  C  =  Q  /  /k*  1  . 

o  o  p  O  r  o 

Functions  of  T 

Parameter  defined  after  equation  22 

Specific  enthalpy  change  during  vaporization 

Electrical  current 

Exponent  of  variation  of  £  with  T 

Thermal  conductivity 

Length  of  conductor  or  conduction  path 

Lengths  of  cold  and  warm  parts  of  two-stage  cooled  conduction  path 

Exponent  of  variation  of  fl  with  T 

Mass  flow  rate  of  vapor 

Number  of  radiation  shields 

Exponent  of  variation  of  k  with  T 

/ - 2 - ' 

Parameter,y  (2  y  /  d  )  -  1,  or  else  difference  n  -  m 
Heat  flux 

Ratio  of  absolute  temperatures,  Tq  /  Tq  . 

Ratio  of  absolute  temperatures,  Tq  /  Tj  . 

Absolute  temperature 
Rate  of  refrigeration  work 

Length  coordinate  of  lead  or  conduction  path  measured  from  the  cold  end 


Vtf 


a,P 

r 


© 

x 

P 

<h 

<P 

r 


Parameters 

Parameter,  =  P  ^  /  2  (kA)ATQ  . 

Emissivity 

Ratio  of  actual  COP  to  Carnot  COP  in  refrigeration  system 

Dimensionless  absolute  temperature,  T/Tq 

Lagrange  multiplier  or  else  ratio,  I  ^  /  Iq 

Electrical  resistivity 

Stefan-Boltzmann  constant 

Rate  of  heat  removed  per  unit  length  at  x  and  T 

Rate  of  heat  removed  per  unit  temperature  change  along  lead 
conduction  path  at  x  and  T 


Subscripts 

a  Ambient  or  condition  at  which  heat  is  rejected 

o  Condition  at  lowest  temperature  in  the  system 

opt.  Optimum  value  or  condition 

1.  Condition  at  intermediate  cooling  temperature 
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Fig. 17.  Refrigeration  work,  W,  for  cooling  a  cryogenic  electrical 
lead  with  resistivity  constant  resistivity,^ ,  and  and  k 

1.  optimum  cooling 

2.  vapor  cooling,  f  -  foptinum 

3.  vapor  cooling,  f  ■  1 

4.  vapor  cooling,  f  ■  0.5 

5.  cooling  only  at  T 
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Fig.  18.  Refrigeration  work,  W,  for  cooling  a  cryogenic  electrical 
lead  with  resistivity  p 

1 .  optimum 

2.  vapor  cooling,  f  ■  1 

3.  cooling  only  at  Tq 


=  >o  (T/T  ),  and  r\  and  k  constant 
f  ft  61 

cooling 


Fig.  19.  Refrigeration  work,  W,  for  cooling  a  cryogenic  conduction  path 

with  conductivity,  k  »  k  (T/T '  )  and  ratio  of  refrigerator  COP  to 

o  ^  a 

Carnot  COP,  =»  n.(T/T.)  • 

6  a 

1 .  optimum  cooling 

2.  vapor  cooling,  f  =  1 

3.  optimum  two  stage  cooling 

4.  cooling  only  at  Tq 


Fig.  20.  Refrigeration  work,  W,  for  two  stage  cooling  at  and  Tq 

with  conductivity,  k  =  k  (T/T  )  and  ratio  of  refrigerator 

^  a  o.3 

COP  to  Carnot  COP,  ri  *■  r]  ( T  / T  )  ’  ,  showing  variations 

&  ft 

with  s. 


1 

! 


,  Fig.  21.  Refrigeration  work,  W,  for  cooling  a  cryogenic  conduction  path 

I  with  conductivity,  k  =  k  (T/T  )  and  constant  tj. 

JL  &  & 

1.  optimum  cooling 

12.  vapor  cooling,  f=  1 

3.  optimum  two  stage  cooling 

4.  cooling  only  at  TQ 

I 


Fig.  22.  Refrigeration  work,  W,  for  two  stage  cooling  at  and  T^, 
with  conductivity,  k  =  k  (T/T  )  and  constant  t),  showing 

a  3. 


variations  with  s. 
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Fig.  23.  Refrigeration  work,  W,  for  cooling  a  cryogenic  conduction  path 
with  conductivity,  k,  and  constant. 

1.  optimum  cooling  (also  vapor  cooling, 

f  _  f  =1) 

optimum 

2.  optimum  two  stage  cooling 

5.  cooling  at  Tq  only 
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I.  INTRODUCTION 


Superconductors  that  carry  very  large  currents  in  high  magnetic  fields  have  been  achieved 
in  the  last  two  years.  Quickly  following  the  finding  by  Kunzler  and  his  associates. 
Reference  1,  that  Nb^Sn  cored  wires  carried  more  than  10^  amperes  per  square  centimeter, 
other  investigators  began  developing  Nb-Zr  wire  and  applying  both  to  magnet  construc¬ 
tion.  Now  we  manufacture  solenoids  that  have  flux  densities  of  well  over  50,000  gauss. 
Reference  2.  What  follows  describes  one  magnet  development  and  sketches  the  engineer¬ 
ing  problems  and  their  solution  for  large,  high  field  magnets.  The  design  is  to  serve  a 
plasma  acceleration  experiment. 


Magnetohydrodynamic  space  propulsion  or  power  generation  systems  requires  electro¬ 
magnets,  which  are  normally  very  heavy  and  require  a  great  deal  of  power.  Suppose 
that  exclusive  of  a  magnet  such  a  system  had  a  weight  W  per  kilowatt  of  power  either 
produced  or  used  in  plasma  acceleration.  The  power  for  a  conventional  magnet  is 
studied  in  Appendix  A. 


Pe 


-0.46 


Eq.  (1) 


And,  the  power  in  say  an  MHD  generator  for  constant  velocity  u  and  conductivity 
is:  Reference  3. 


Eq.  (2) 


When  Pe  has  to  be  supplied  by  output  P  and  the  weight  of  the  magnet  is  added  to  the 
system  weight  WP,  the  weight  power  ratio  can  be  found  making  use  of  Equations  (1) 
and  (2)  to  eliminate  dimensions: 
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Eq.  (3) 


Weight 
Net  Power 


W  + 


-  - 2~5 - 

4ir^pe(pgAh)(V°2)  0,46 

2 

H  XaUr|e0-r|e)P 
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Typical  values  might  be  a,  10  mhos/m;  B,  50,000  gauss;  u,  10  rr/sec;  ^  ,  0.  8;  p  , 

-8  4  3  6 

2  x  10  ohm  m;  (p  Ah),  4  x  10  joule^/m  ;  X,  1;  and  say  W  of  10  ll/kw,  and  P, 

9 

12,  000  kw.  Then: 


Weight 
Net  Power 


_  10  +  0.  25(V°2) 

i  /A/  2. -d.  46 

1  -  (A/a  ) 


(ll/kw) 


Eq.  (4) 
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Optimum  A/a  is  near  10  so  20%  of  total  weight  is  in  the  magnet  and  35%  of  P  goes 
to  the  magnet.  A  very  similar  consideration  applies  to  plasma  propulsion.  This  very 
simplified  estimate  ignores  magnet  joule  heat  rejection.  At  P  equal  1  megawatt,  no 
net  power  could  be  obtained. 


But,  with  the  use  of  superconductors,  high  j  leads  to  reduced  conductor  weight;  no  need 
for  iron;  elimination  of  Pe  except  for  a  modest  amount  of  refrigeration  power;  and  an 
ability  to  develop  strong  fields  which  shorten  length  |  to  reduce  viscous  and  heat  transfer 
irreversibilities  or  decrease  electrical  plasma  dissipation  fraction  1  ~r)e. 


II.  SUPERCONDUCTOR  PROPERTIES 


Although  many  metals  are  superconductors,  at  temperatures  below  20  K,  recent  attention 
has  been  focused  on  such  intermetal  lie  compounds  typified  by  Nb^Sn,  V^Si  and  V^Ga 
and  interstitial  alloys  typified  by  body-centered  cubic  systems  around  Nb  in  the  periodic 
table,  especially  Nb-Zr.  Unlike  ideal  superconductors  which  have  critical  fields  of  the 
order  of  10,000  gauss  and  less,  these  materials  exhibit  mixed  state  behavior  where  flux 
does  begin  to  penetrate  wire  portions  that  become  normal  conductors  at  fields  above  the 
thermodynamic  critical  field,  but  zero  electrical  resistance  persists  to  a  much  higher  field. 
Critical  current  density  and  fields  are  shown  in  Figure  1  which  is  taken  from  data  given  by 
Kunzler,  Reference  4,  for  Nb^Sn  and  Hulm,  Reference  5,  for  Nb-25%  Zr.  Magnets  up 

to  200, 000  gauss  appear  to  be  possible  using  Nb.Sn  and  up  to  90, 000  gauss  using  Nb-Zr. 

J  r  T  n2 

Variation  of  critical  field  with  temperature  follows  approximately  =  Hq  [1  -  ( ^  )| 
where  critical  temperature  Tc  for  Nb^Sn  is  18°K  and  11°K  for  Nb-25%  Zr.  The  curves 
of  Figure  1  shift  up  and  right  for  lower  temperature  than  4.  2°K;  down  and  left  for  higher 
temperature.  Wire  performance  in  coils  using  maximum  fields  on  the  center  inner  wind¬ 
ing  of  solenoids  for  data,  shows  reduced  critical  currents  from  those  in  short  samples  in 
all  experiments  to  date.  This  degradation  in  Nb-Zr  is  least  for  50%  Zr,  Reference  5, 
but  short  sample  currents  are  also  less  than  for  lower  Zr  content.  Optimum  j  occurs  with 
25  to  33%  Zr.  The  degradation  appears  to  result  from  a  proximity  effect  of  the  near 
range  individual  conductor  fields  on  each  other  and  the  resultant  current  distribution  in 
wire.  Current  density  also  increases  with  cold  work  reduction  in  wire  making.  Because 
it  is  ductile  we  have  used  Nb-25%  Zr  wire,  which  has  been  hard  drawn  from  0.50  inch 
diameter  rod  to  0.010  inch  diameter  wire.  Niobium-tin  wire  is  brittle,  making  coil 
fabrication  difficult.  While  the  metallurgical  structure  gives  the  mixed  state  supercon¬ 
ductor  properties,  it  aiso  gives  high  normal  state  superconductor  resistivities  which  lead 
to  one  of  the  major  engineering  problems,  coil  protection. 
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M  (KILOGAUSS) 


Figure  1  -  Critical  Current  Density  in  0.006  in.  Nb„  with  0.011 

in.  Diameter  Nb  Sheath  and  .  010  in.  Diameter  Nb- 25% 
Zr  Wire  at  4.  2°K 
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III.  MAGNET  DESIGN 


MAGNET  SHAPE  AND  FIELD 

Magnets  for  application  to  linear  MHD  devices,  must  produce  somewhat  uniform  cross 
fields  in  a  volume  which  is  roughly  square  and,  to  avoid  irreversible  end  effects,  with 
a  length  to  width  ratio  greater  than  unity.  The  coil  shape  shown  in  Figure  2  was  used 
because  it  minimizes  weight,  can  produce  useful  fields  near  the  superconductor  critical 
field,  has  relatively  good  uniformity  of  field,  and  it  can  be  and  has  been  fabricated. 

The  formulas  in  Appendix  A,  for  case  5,  giveit  Bq/ pjaX  equal  1. 56  for  a  of  2.  2, 

of  1. 1  and  (3  for  a  midplane  split  equal  0.  05.  This  infinite  conductor  result  can  be 

increased  by  7%  to  account  for  finite  length  and  end  conductors.  The  correction  was 

calculated  for  a  similar  shape  using  right  angle  bends  and  rectangular  sections.  Then 

8  2 

with  X  equal  0.  5,  a  equal  2.  2  inches  and  j  of  2. 76  x  10  amp/ m  for  14  amps  in  0. 010 

2 

inch  wire,  a  field  of  5. 14  webers/ m  would  be  developed. 

If  a  simple,  split  solenoid  pair  whose  shape  has  been  optimized  were  used  to  obtain 
the  same  field  over  an  effective  length  of  10  inches,  its  inside  diameter  would  have  to 
be  7. 58  inches  and  2.  3  times  the  weight  of  wire  would  be  used,  because  of  the  large 
split.  Greater  duct  length  to  width  ratios  would  increase  the  weight  ratio  even  more. 
Worse  yet  would  be  field  non-uniformity  which  for  example  would  also  cause  a  maximum 
field  on  superconductors  to  be  1. 37  times  center  field,  limiting  MHD  fields  to  73%  of 
critical  field. 

Coil  winding  techniques  were  investigated  and  one  selected  which  after  development 
produced  the  shape  in  Figure  2,  which  is  essentially  like  a  race  track  rolled  on  to  a 
cylinder. 
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Figure  2.  Superconducting  Magnet  Winding  and  Coil  Support 


Magnetic  energy,  T  or-Tj  Li  ,  for  the  design  has  been  calculated  from  Equation  (B-2) 
in  Appendix  B  for  equivalent  rectangular  conductor  bundles  (a  =  2,  p  =  1).  In  view 
of  the  work  by  Arnold,  Reference  6,  the  result  obtained  using  a  10  inch  effective 
length  should  be  quite  accurate.  There  are  about  87,000  joules  when  Bq  is  5 
webers/m^  and  therefore,  L  is  940  henries.  Dissipating  this  large  amount  of  energy  is 
accompanied  by  a  magnet  protection  problem  should  the  conductors  enter  the  normal 
state. 


COIL  PROTECTION 


Though  it  is  possible  to  operate  superconducting  magnet  systems  reliably,  but  by  some 
accident,  e.g.  loss  of  refrigeration,  the  conductors  enter  the  normal  state,  then  in 
spite  of  dissipation  of  large  amounts  of  magnetic  energy,  coils  must  be  protected  from 
destruction  of  the  continuity  of  the  superconductor.  A  simply  wound  coil  having  a  por¬ 
tion  resistive  has  the  circuit  shown  in  Figure  3(a).  If  externally  supplied,  by  a  battery 
and  control  resistor,  still  these  are  not  significant  circuit  elements  in  a  transient.  Ther¬ 
mal  energy  increases  temperature  in  normal  portions  and  diffuses  into  superconducting 
regions  causing  propagation  and  growth  of  R  which  is  proportional  to  the  size  of  the 
normal  region  and,  as  resistivity  changes  slightly,  its  average  temperature.  Discussion 
of  thermal  propagation  appears  in  Appendix  C.  For  constant  resistivity  and  turn  induc¬ 
tances  the  transient  is  summarized  below: 

L  4-  +  f(T)R  i  =  0,  i(o)  =  i  ,  f(T  >  t  )  =  1  Eq.  (5) 

dT  o  o  o 


-R 


i 


=  e 


‘f  dT 

°  ,  0<T  <T 


\  Eq.  (6) 


f  i  f  (t  “  To> +  f  fd*rl 


-r-  =  e  T  l 


/  T  >  T 
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NORMAL  SUPERCONDUCTING 


UNPROTECTED  COIL  LOW  RESISTANCE 


(a)  (b) 

unouAi  ci  iDCDmuniv'Tiun 


COIL  WITH  MANY  SHUNTED  SECTIONS 

(d) 


Figure  3  -  Superconducting  Coil  Circuits  During  Normal  Transients 
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Eq.  (7) 


2V  r  R 
n  _  o  0 
TFT'MTT 


2i  R 
o  o 

“TT" 


0  -Oe 


0  <  T  <  T 

0 

Largest  voltage  gradients  occur  between  ends  of  pairs  of  the  N  layers  wound  back 
and  forth  in  a  coil.  That  layer  to  layer  voltage  for  the  normal  region  is  given  by 
Equation  (7)  which  shows  that  it  is  largest  when  two  whole  layers  are  normal,  and  the 
fraction  f  is  small  so  that  voltage  is  twice  i  times  layer  normal  resistance.  Inductive 
gradient  in  the  superconducting  regions  must  be  smaller  since  the  factor  1-f  is  re¬ 
placed  by  f  in  Equation  (7)  and  when  f  is  large,  the  exponential  is  small.  An  entirely 
normal  coil  has  no  voltage  gradients  by  this  simple  theory.  Should  insulation  break¬ 
down  and  arcing  occur,  local  heating  anneals  wire  and  destroys  or  weakens  super¬ 
conducting  properties.  Maximum  voltage  2  i^R^/N  can  be  alternately  written  as 
2jpe^  (layer  volume)/  (wire  cross-section  area).  Clearly  the  protection  problem 
arises  from  extremely  high  supercurrent  densities,  rather  high  normal  resistivity  and 
small  wire  diameter,  for  practical  layer  volumes.  The  MHD  magnet  design  would  have 
a  possible  18,000  volts  layer  end  to  layer  end  were  no  protective  techniques  used. 


A  way  that  coil  current  and  therefore  high  voltage  may  be  reduced  is  by  coupling 
magnetically  from  a  short  time  constant  (L/fR^)  normal  superconductor  system  to  a 
secondary  which  has  a  significantly  longer  time  constant  (L^l^)  because  it  has  the 
low  resistivity  of  copper  or  aluminum  at  cryogenic  temperatures.  The  secondary  of 
Figure  3(b)  now  alters  the  problem  of  Equation  (5)  to: 

cT  di  ' 

L  +  f(T)RQi  +  Mg|-  =  0,  i  (o)  =  io 

di«  f 


+  i2R2  =  *2^  =  0 


Eq.  (8) 
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Suppose  that  this  system  is  closely  coupled  (k  =  M  /  L  Lj  =  1)  and  that  as  an  approxi- 

'  (8) 

t/t, 


motion  f  (t)  =  t/  t  ,  then  solution  of  Equation  (8)  gives: 


2V  2R 


2i  R  (I  -  t/t  Je 
n_  o,  oo  _ o_ 

JFT~FT  0  0  1  — FT" 


l  +  -i  i 


OJ 


'-Vo' 


Eq.  (?) 


T  <  T 


Clearly  from  Equation  (9)  for  a  secondary  to  be  effective  t  ^  should  be  much  greater  than 

t,  and  Vt,t  .  Then  voltages  are  reduced  quickly  in  a  characteristics  time  t,t  /t« 
1  1  o  '  1  o  2 

T  o 

by  a  factor  of  1  +  f  — -  .  However  if  there  is  not  close  coupling,  a  solution  of  Equation 

T1 

(8)  by  a  series  in  time  gives: 


2V  2i  R 
n  _  o  o 

fN  “ TT 


1  - 


2t.t  (1  -  k^) 
i  o 


+  . , 


Eq.  00) 


When  T2  is  truly  very  large  so  that  additional  terms  in  Equation  (10)  are  unimportant, 
then  voltages  may  not  be  reduced  nearly  instantaneously  in  time  t^t  /  t ^  but  decay 
/  2 

in"\/2T,T  (1  — k  ).  In  either  case,  final  slow  decay  occurs  with  a  t«  time  constant. 

v  I  o  z 

Similarly  when  a  fraction  f  of  a  coil  goes  normal  suddenly,  a  reduction  of  Equation  (7) 

°  2 
by  1  +  f  t 2/  t ^  occurs  instantaneously  for  close  coupling.  When  k  is  not  one, 


reduction  occurs  with  a  time  constant  (1  -  k  )  T,/f  . 

1  o 


Another  means  of  protection  is  resistance  shunting  of  several  coil  sections.  A  coil  with 
just  two  shunted  sections  where  one  is  normal  is  shown  in  Figure  3(c).  The  governing 
equations  are: 


h  -3T  *  m7T  +  <MVi  +  Rs(i,  -  V  ‘  ' 

di!  di2  ) 
M7T  +L2lT  +  Rs(i2‘  'l)_0'  *2  (o)  -  ’o 


Eq.  01) 
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For  close  coupling  and  f  (t)  =  t/t  /  the  solution  is: 

o 


2V  2i  R  r 

n  o  o 

,  t/to  1 

-t/t2 

e 

7n - FT-  [ 

_ i 

> 

+ 

u- 

1 

h-  -1 

T,  T  1 

l  1  0  J 

1  TlT°l 
r? 
t2 

'  T<flT0 


Eq.  12 


A  comparison  of  Equations  (1  2)  and  (9)  shows  that  shunted  behavior  is  very  nearly  like  that 
for  normal  secondaries  where  now  the  secondary  is  composed  of  and  R^.  Since  the  time 

constant  t9  may  be  made  arbitrarily  large,  the  characteristic  time  for  current  and  voltage 

1  _ - 2 

decay  to  the  slow  transient  is  important.  Here  the  time  constant  is  2t«t  (1  -  k  )  for 
2  \  oo 

f  =  t/tq  and  (1  -  k  )/f^  for  instantaneously  appearing  resistance  fQR  .  The  last 

results  were  obtained  from  series  solution  of  Equation  (12). 


The  usefulness  of  shunting  sections  of  a  coil  is  the  rapid  current  decay  in  sections  hav¬ 
ing  normal  wire,  accompanied  by  rapid  current  increase,  above  iQ,  in  the  most  closely 
coupled  superconducting  sections.  Analysis  of  the  voltage  across  the  terminals  of  the 
"k"  section  of  Figure  3(d)  shows: 
n 

I  k:  dl:  /  i“  1  1  si  .  \  R 


V  di.  i  £  !  'A 

r-i  '^  +  5kRk= w* -(■!?“ 


r*  R 

il  si 

i  =  1  1 


L-M  "* 


Eq.  (13) 


In  the  limit,  as  n  approaches  an  infinitely  large  number  Equation  (13)  becomes: 


J  M(k, 


(X'T>dX-'° 
'  3t  n 


[  l(X,T)R(X)dX 

- - -  +  i  (k , t  )  [R (k)  +  R$(k)j  =  0  Eq.  (14) 


/ 


R$  (X)dX 


Since  the  integrals  are  continuous,  while  R(j)  may  be  discontinuous  at  a  division  of 
normal  and  superconducting  sections,  for  example  between  the  k  and  k+  1  sections,  the 
current  ratio  from  Equation  (14)  must  be: 
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Eq.  05) 


*(k  +  1). 
'(k) 


R(k)  +  R  (k) 
R$(k  +  1) 


For  low  resistance  shunts  the  ratio  may  be  very  large  except  that  it  may  be  assumed 
that  section  k  +  1  must  itself  become  normal  for  currents  greater  than  i^  (critical 
current).  The  limiting  process  leading  to  Equations  (14)  and  (15)  implies  close  coupling 
as  for  a  finite  size  solenoid  with  infinitesimal  sections,  which  cannot  be  realized  in 
practice.  But  if  in  the  system  of  Figure  3(c)  there  is  not  close  coupling,  series  solution 
of  Equation  (11)  gives: 


i0  f,  R  k 

J  =  i  .  a- 


Tjr^n-k2) 


T  + 


Eq.  (16) 


From  Equation  (16)  where  f,  R  appears  instantaneously  because  current  exceeds  i 

•  •  vV^i-n 


then  the  second  section  becomes  normal  in  time  ( 


i  - 

c 


In  an  n 


Vv 

section  coil  where  all  sections  are  the  same  and  coupling  effects  are  limited  to  one 
normal  and  one  superconducting  pair  (where  k  is  largest)  a  representative  time  for  going 

2  .  What  has  been  described  amounts  to  a  fast,  non¬ 


normal  is:  i  i  — i 
c  O; 


I 


Tll 


1  -k 

F 


thermal  propagation  of  a  normal  front.  Behind  the  wave,  sections  have  dissipated 
energy  in  the  characteristic  time  found  from  Equation  (16).  Temperature  rise  is  there¬ 
fore  of  the  order  of: 


AT;s  'o^'c  V 


s  L1L2  (1  -  k2) 


CM, 


Eq.  07) 


It  is  conceivable  that  when  i^  decays  sufficiently,  the  section  may  still  be  at  a  low 
enough  temperature  that  it  may  become  superconducting  again.  However  it  would  be 
just  at  the  critical  condition  and  any  magnetic  coupling  would  increase  current  and 
return  it  to  the  normal  state. 


The  coil  described  in  this  paper  is  protected  by  96  shunts  and  46  interleaved  copper 
secondaries.  Were  the  coil  not  shunted,  even  though  rather  close  coupling  with 
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secondaries  would  be  achieved  (k  =  0.975)  and  with  copper  occupying  10%  as  much 

space  as  superconductor,  the  reduction  in  maximum  voltage  shown  in  Equation  (9) 

would  not  be  significant  for  very  small  f  =  t/ t  .  However,  by  subdividing  the  coil 

with  shunts  so  that  one  secondary  couples  effectively  with  one  normal  section  of  96, 

reduction  may  be  attained.  The  shunting  permits  the  current  of  one  section  to  couple 

independently  with  any  normal  secondary  or  superconducting  section.  For  two  layers 

wholly  normal  in  a  section  and  with  t  ,  /t  ,  =  20  for  copper  at 

'  secondary  normal  section  rr 

20  K,  at  most  2500  volts  layer  to  layer  could  occur  even  without  section  to  section 
coupling.  The  way  in  which  thermal  propagation  of  a  normal  front  occurs  in  a  section 
is  complicated  by  cooperative  heating  in  secondaries.  A  section's  tq  is  conservatively 
estimated  as  10  ^  seconds  (see  Appendix  C).  Layer  to  layer  voltage  becomes  a  problem 
only  when  t/tq  becomes  nearly  two  divided  by  the  layers  in  a  section.  Then  if  there 
is  not  close  coupling  the  time  constant  of  Equation  (10),  for  this  design,  need  only  be 
like  1/4  t  *  Since  is  much  less  than  tq,  coupling  need  not  be  good.  Then  the 
reduction  predicted  by  Equation  (9)  is  attainable  before  large  voltages  appear.  When 
two  layers  become  normal  suddenly,  there  is  the  time  (l-k^)  Tj/f  (equal  8x10^ 
seconds)  for  which  reduction  of  voltage  does  not  occur.  However,  capacitive  energy 
storage  in  the  dielectric  between  layers,  prevents  voltage  rise  for  of  the  order  of  10  ^ 
seconds,  so  protection  is  achieved  under  this  severe  condition  too. 

The  shunt  action,  provided  adjacent  sections  could  remain  superconducting,  would 
produce  a  250  times  reduction  of  maximum  layer  to  layer  voltage  by  Equation  (12)  or 
(15).  Since  in  the  design  for  maximum  field  they  cannot  remain  superconducting,  the 
useful  effect  of  shunts,  by  themselves,  is  to  propagate  the  normal  region  in  .19  x  10  ^ 
seconds  per  average  section.  In  the  design,  shunted  section  time  constants  initially 
are  longer  than  secondary  time  constants  so  that  the  propagation  is  not  hindered  by 
copper  secondaries,  at  least  for  comparable  coupling.  Temperature  rise  in  this  time 
is  of  the  order  of  10°K  by  Equation  (17),  insuring  that  sections  remain  normal.  If  each 
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section  was  made  completely  normal  by  current  increase,  no  voltage  problem  would 
exist  and  copper  secondaries  would  be  redundant  except  to  slow  final  field  decay. 
However,  the  secondaries  do  provide  protection  should  only  a  fraction  of  each  coil 
become  normal,  and  the  preceding  discussion  of  secondaries  again  applies. 

Another  means  of  protection  which  has  not  been  used  is  that  where  the  superconductor 
is  clad  or  plated  with  copper  and,  therefore,  has  a  much  reduced  voltage  problem  by 
Equation  (7).  In  effect,  this  copper  is  equivalent  to  a  system  of  shunts  for  current 
continuity  and  secondaries  for  magnetic  energy  conservation.  However,  one  further 
problem  does  exist.  It  is  the  final  distribution  of  thermal  energy  in  a  coil.  Were 
thermal  propagation  of  normalcy  slow,  high  temperatures  would  occur  in  the  first 
superconductor  to  become  normal,  especially  in  large  coils.  The  shunt  system  over¬ 
comes  this  by  the  extreme  rapidity  with  which  it  spreads  the  normalcy.  Then  only  in 
extremely  high  energy,  large  magnets  does  the  resulting  uniform  distribution  of  thermal 
energy  cause  temperatures  high  enough  to  destroy  the  superconductor  characteristic. 

While  the  foregoing  contains  many  approximations,  such  as  the  linear  thermal  propa¬ 
gation  or  sudden  transition  approximation  and  the  assumption  of  constant  properties, 
the  conclusions  and  descriptions  remain  valid.  A  computer  program  including  non¬ 
linear  thermal  transfer,  changing  critical  current  with  field  and  temperature,  exact 
coupling  coefficients  and  circuit  analysis  yields  us  detailed  behavior  and  exact 
solutions. 

SUPPORT  OF  MAGNETIC  FORCES 

In  the  system  of  conductors  shown  in  Figure  2,  the  net  result  of_[  x  B  body  forces,  is  to 
tend  to  enlarge  the  coil,  to  separate  the  sides  by  a  force  of  68,500  lb.  found  by 
Equation  (D-2  of  Appendix  D),  where  a  =  2,  J3  =  1 ,  and  to  separate  ends  by  a  force  of 
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77,000  pounds.  Average  stress  on  the  coil  form  sides  is  quite  accurately  found  from 
Appendix  D,  Equation  (D-l)  to  be  1360  psi  and  varies  less  than  17%  over  the  height 
of  the  side.  The  maximum  compressive  stress  between  halves  is  only  a  little  more, 

1640  psi.  For  still  larger  sizes  and  fields,  stress  increases.  With  the  fully  potted  coil 
structure  used,  epoxy  insulation  will  not  be  damaged  by  these  stresses. 

The  magnet  described  here  is  uniike  a  simple  solenoid  where  the  coil  supports  itself 
with  a  stable  conductor  bundle  pinch  and  hoop  stress  jrB  (for  a  free  wire  turn  of  radius, 
r,  with  an  axial  field  B  at  the  wire)  constraining  coil  enlargement.  Rather  the  extra 
constraining  coil  form  structure  shown  in  Figure  2  is  used. 

e 

Plate  loading  is  reduced  using  end  and  side  ribs  which  also  act  as  beams  to  communicate 
load  to  axial  and  cross  tension  bolts.  Side  and  end  flanges  are  bolted  to  connect  the 
beams.  The  side  flanges  also  act  as  beams  to  carry  side  loading  to  circular  end  plates 
which  acts  as  connecting  tension  members.  End  loading  is  also  carried  by  the  side  coil 
form  and  by  the  flanges  in  superimposed  tension.  All  members  are  stainless  steel  designed 
for  50,000  psi  maximum  stress  at  4.  2°K,  except  where  there  is  redundant  sharing  of  load. 
For  example  both  side  flanges  plus  coil  form  and  longitudinal  axial  bolts  can  support  the 
end  loading  by  themselves. 

ELECTRICAL  CONNECTIONS 

Over  200, 000  feet  of  wire  was  used  in  the  coils.  Although  up  to  18, 000  foot  lengths 
of  Westinghouse  produced  wire  were  supplied,  the  average  single  conductor  length 
after  coil  fabrication  was  much  less.  Superconducting  junctions  were  made  by  pressing 
twisted  wire  ends  together  between  the  bottom  of  a  milled  slot  in  a  stainless  steel  bolt 
and  a  foot  tightened  by  a  nut  on  the  bolt.  These  junctions  were  tested  up  to  20, 000 


gauss  fields  and  found  to  have  nearly  the  same  critical  current  as  the  wire.  Junctions  are 
clamped  over  insulated  fingers  on  a  copper  conduction  cooling  strip.  In  order  to  divide 
the  coil  into  eight  sections  for  external  supply  of  current  into  individual  sections  and 
also  to  furnish  different  currents  to  the  sections,  seven  three  wire  superconducting  joints 
are  used.  The  third  wire,  a  superconducting  lead,  is  carried  through  vacuum  seals  into 
wells  on  liquid  helium  tanks  where  connection  is  made  to  copper  wire.  The  nine  current 
joints  are  simply  copper  blocks  clamping  the  superconductor. 

Shunts  are  conventional  2  watt  wire  wound  resistors,  which,  because  resistance  wire 
transmits  heat  well  to  high  specific  heat  plastic,  can  absorb  100  joules  of  energy  without 
failure.  During  a  normal  transient  they  can  carry  14  amperes  in  10  ohms  for  .05  sec.,  a 
time  that  is  longer  than  that  expected  for  shunt  currents  to  flow.  They  connect  at 
junction  blocks  where  either  superconducting  joints  or  continuous  conductors  are  brought 
out  of  the  coil. 

Figures  5  and  6  show  the  wire  carrier  used  and  the  superconducting  junctions. 
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IV.  CRYOGENIC  SYSTEM 


GENERAL  CONSIDERATIONS 

Superconductors  need  to  be  maintained  at  the  normal  saturation  temperature  of  helium, 

4. 2°K.  Lower  temperatures  increase  refrigeration  requirements  markedly  without 
allowing  much  increase  in. critical  currents.  Higher  temperatures,  up  to  a  maximum  of 
11°K,  resultin  progressively  smaller  current  densities.  Refrigeration  work  to  remove 
Q  units  of  heatct  temperature  T  and  reject  at  Tr  is: 

W  =-  Q  Eq.  (18) 

*i  1 

Here  ^  is  the  ratio  of  the  coefficient  of  performance  obtained  to  that  of  a  Carnot  cycle. 
For  example  for  a  helium  cycle  ^  is  between  7  and  20%  while  rejecting  at  an  Earth 
ambient  of  300°K.  For  the  lower  value,  Q/W  is  0. 001 ;  that  is,  1000  watts  of  power 
are  required  to  remove  a  1  watt  heat  flux  into  the  cold  region.  The  need  for  minimizing 
heat  transfer,  especially  to  the  lower  temperatures  is  well  demonstrated  by  this  con¬ 
sideration  whether  it  be  for  the  experimental  development  described  here  or  for  magnet 
cryogenics  in  space. 

Heat  transfer  occurs  by  radiation  from  warm  walls  to  the  magnet,  by  conduction  in 
residual  gas  unless  vacuum  pressures  less  than  10  ^  mm  mercury  are  used  to  make  this 
heat  flux  negligible,  by  conduction  in  magnet  mass  and  reaction  supports  and  conduction 
and  joule  heating  in  electrical  current  ond  instrumentation  lead  conductors.  Whether  the 
object  is  to  minimize  weight  and  power  of  a  space  system  refrigerator  or  to  conserve 
liquid  helium  in  experiments,  refrigeration  to  remove  as  much  heat  as  possible  at  temper¬ 
atures  above  4. 2°K  has  the  advantages  implied  in  Equation  (18).  Since  thermal  conduc¬ 
tivity  in  supports  and  radiation  fluxes  decrease  rapidly  with  temperature,  staged 
refrigeration  is  especially  practical.  Although  it  does  not  result  in  an  optimum  temper¬ 
ature  distribution  of  refrigeration,  vapor  cooling,  using  enthalpy  changes  in  the  vapor 
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as  well  as  evaporation  is  a  good  and  simple  technique,  it  may  be  applied  to  both  con¬ 
duction  and  radiation  transfers  through  vapor  cooling  of  supports,  radiation  shields, 
electrical  leads,  etc.  For  example,  for  constant^  and  properties  k  and  c^.  Equation 
(E-5)  of  Appendix  E  together  with  a  liquification  process  following  Equation  (18)  gives 
refrigeration  power  for  vapor  cooling  a  heat  path  os: 


P  Vapor  .1  kA  , 

c  (T  -T) 

1  +  P,  ° 

c  rVa-T 

,T  T  l' 

dT  +  hf 

r  Cooling  rj  c^x^  n 

h. 

L  fg  J 

LpJ  (  T 

«»Pp 

Without  vapor  cooling  the  result  would  be: 


Eq.  (19) 


Pr  No  =  1  kA  (Ta-Ts)  Eq.  (20) 

Vapor  q  x  T 
Cooling  °  5 


Refrigeration  power  without  vapor  cooling  would  be  four  times  that  with  vapor  cooling 
considering  transfer  between  4.  2°K  and  300°K  with  helium.  A  similar  consideration  in 
the  case  of  cooling  electrical  leads  yields  vapor  cooling  improvement  by  a  factor  of 
about  three  in  a  comparison  of  optimum  cases. 


THE  MAGNET  DEWAR 

Figure  4  shows  a  schematic  of  the  dewar  design.  The  vacuum  vessel  is  stainless  steel 
made  in  two  sections.  A  3.  9  inch  I.  D.  by  18.  6  inch  long  center  tube  about  the  MHD 
duct  space,  a  head  end  plate  and  a  lead-in  neck  comprise  one  part.  The  other  is  a 
20.5  inch  diameter  outer  cylinder  and  rear  end  plate  which  has  "O"  ring  seals  at  a  head 
end  flange  and  at  the  rear  end  of  the  center  tube.  This  arrangement  allows  relatively 
easy  assembly. 


A  second  plate  at  the  rear  end  and  a  thin  stainless  steel  cylinder  form  a  part  of  the 
support  structure,  acting  also  as  intermediate  radiation  shields  and  low  emissivity  surfaces. 
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NITROGEN  TANKS 


JRE  VESSEL  TO  NITROGEN 
SUPPORT  CABLES 


ivAi/i  r  nnrJ  ^ncnAnnion  System 


Four  0. 10  inch  diameter  stainless  steel  cables  at  each  end  support  and  center  the 
nitrogen  system.  Eight  additional  cables  prevent  axial  displacement. 

The  nitrogen  system  has  two  copper  end  radiation  shields  with  heavy  outer  rings  to  which 
the  support  cables  are  fastened.  A  stainless  steel  inner  cylindrical  shield  is  plated  with 
copper  strips  arranged  to  conduct  heat  to  the  end  shields  but  to  offer  a  high  resistance 
path  for  induced  eddy  currents  which  may  deform  it.  An  outer  shield  is  a  thin  copper 
sheet.  Heat  is  conducted  to  three  interconnected  cylindrical  tanks  at  the  bottom  which 
hold  2.9  liters  of  liquid  nitrogen.  A  single  0.010  inch  wall  stainless  steel  fill  and  vent 
tube  connects  with  the  top  of  the  neck.  Inside  it  is  a  small  teflon  tube  used  for  filling. 
Eight  adjustable  support  cables  run  diagonally  from  the  nitrogen  temperature  ring  to  the 
coil  form  flanges. 

Liquid  helium  is  stored  in  two  tanks  holding  3.65  liters  which  are  mounted  on  top  of  the 
coil  form.  Stainless  steel  sheet,  plated  with  copper  in  strips,  surrounds  the  magnet. 

They  conduct  radiated  heat  to  the  helium  reservoirs  at  right  angles  to  electric  fields 
which  could  induce  eddy  currents  during  magnet  normal  transients.  These  shields  are 
needed  since  many  of  the  magnet  surfaces  offer  high  emissivities.  Thin  wall  tubes  run 
parallel  from  the  neck  to  the  head  end  nitrogen  ring  where  they  are  heat  shunted  to 
nitrogen  temperature.  From  there,  one  serves  as  a  fill  line  which  connects  with  a  vapor 
cooling  heat  exchanger,  while  the  other  vents  the  two  tanks  at  the  rear  end. 

Helium  and  nitrogen  temperature  shield  surfaces  were  gold  plated  to  minimize  oxidation 
and  deterioration  of  emissivities.  Nitrogen  shielding  of  radiation  and  conduction  paths 
effects  the  staging  discussed  above.  Vapor  cooling  of  electrical  leads  and  tubing  was 
calculated  as  optimum  using  equation  (E-3  of  Appendix  E)  with  appropriate  convection 
heat  transfer  theory.  In  this  connection,  it  was  noted  that  the  effective  h  was  increased 
by  using  heavily  insulated  lead-in  wires  since  teflon  thermal  conductivity  is  not  so  low 
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that  the  additional  resistance  to  heat  flow  by  insulation  offsets  the  gain  in  convection 
transfer  area.  Nine  current  leads  are  provided  to  the  magnet  and  since  all  would  not 
be  carrying  the  planned  current,  their  length  to  area  ratio  was  chosen  to  give  smallest 
heat  input  to  helium  liquid  at  either  full  or  no  current,  i.e.  xq/A  in  Equation  (E-4  of 
Appendix  E)  large  rather  than  optimum  by  Equation  (E-6).  Table  I  summarizes  the  cal¬ 
culated  heat  transfer  in  the  dewar  for  the  superconducting  magnet. 

MAGNET  COOLING  AND  COOL  DOWN 

The  winding  is  kept  cold  by  conduction  from  insulated  tabs  on  pure  copper  secondaries 
to  boiling  helium  in  blind  tube  wells  communicating  with  the  helium  tanks.  Joint  and 
superconductor  lead  carriers  are  similarly  kept  cold.  A  helium  fill  line  of  stainless 
steel  connects  to  two  vapor  heat  exchangers,  also  cooling  secondaries,  so  that  efficient 
vapor  cool-down  from  nitrogen  temperature  is  effected  as  described  in  Appendix  F.  The 
200  pound  magnet  assembly  would  require  174  liters  of  boiling  helium.  However,  using 
Equation  (F-l)  and  (F— 4)  gives  a  calculated  requirement  of  14.  2  liters  for  this  design, 
exclusive  of  transfer  losses.  The  boiling  tubes  and  vapor  heat  exchanger  may  be  seen 
in  Figures  5  and  6. 

During  electrical  charging,  voltage  for  L  di/dT  also  causes  dissipation  of  heat  in 
secondaries  and  shunts.  In  shunts  for  example,  total  energy  to  charge  to  i  in  time  t 
at  a  constant  rate  is  ,2L  W1  .  .  2V  A  long  charging  time  is  used  to  minimize 

VTT'  1  'o  ’ 

s 

dissipation.  Heat  transfer  calculations  have  been  used  to  determine  that  the  design  will 
not  warm  superconductors  appreciably  during  charging.  Also,  shunts  have  been  chosen 
so  that  sections  have  nearly  the  same  currents  during  charging. 
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Table  I  -  Heat  Transfer  in  Dewar 


I 

I 


Figure  5  - 


ing  Magnet  Bottom  View 


23 


Figure  6  -  Superconducting  Magnet  Assembly 
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V.  INSTRUMENTATION 


In  addition  to  use  of  the  nine  current  leads,  to  the  magnet,  five  voltage  leads  are  pro¬ 
vided  to  monitor  transient  section  voltages.  Copper-constantan  thermocouples  are 
installed  to  measure  nitrogen  tank,  magnet  support  and  helium  vapor  exhaust  temperature. 
Monitoring  the  latter  two  enables  setting  charging  rates  to  obtain  efficient  vapor  cool 
down.  A  nitrogen  level  gauge  and  helium  level  gauge  are  used  to  give  continuous 
readings  of  tank  fluid  level.  These  level  gauges  are  described  in  Appendix  G. 


APPENDIX  A 
Fields 


The  flux  density  normal  to  the  axis  of  symmetry  between  two  infinite  parallel  conductor 
bundles  is  twice  the  sum  of  contributions  by  current  elements  in  one.  (See  Figure  7). 


B 

o 


=  J* 

ir 


Eq.  (A— 1 ) 


The  electrical  power  consumed  by  the  pair  in  resistive  magnets  is: 


Pe  =  2peL 


dA 


Eq.  (A-2) 


A  current  distribution  that  minimizes  Pe  is  found  by  a  variational  method  where,  given 
Bq,  Equation  (A— 1 )  is  then  a  constraint. 


•  —  •  /  X  \ 

'opt.  °  <-7— T> 


x  +  y 


Pe  =  2p  I  (irB  /p)' 


IL  ^ 


dA 


r  rr 


rr 


a 

x  +  y 


4  -dA 


Eq.  (A- 3) 


Eq.  (A -4) 


When  j  is  constant  or  optimum  from  Equation  (A-3)  inspection  of  the  alternative  form  of 
Pe  in  Equation  (A— 4)  shows  that  only  the  integral  in  the  denominator  of  Equation  (A-4), 
which  is  irB^/p,  may  be  maximized  to  minimize  Pe,  as  long  as  area.  A,  is  considered 
to  be  fixed  or  constrained.  The  variational  procedure  maximizing  irB^/ p,  with  A  con¬ 
stant,  gives  an  optimum  limit  in  the  integration  on  y.  The  best  shape  conductor  bundle 
both  for  optimum  j  distribution  and  for  constant  j,  is  a  circle  touching  the  axis.  A 
simpler  shape  is  a  rectangular  one  where  again  a  variational  method  yields  the  optimum 
relation  of  dimensions.  Note  that  the  above  procedure  finds  minimum  power  for  area 
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ELEMENT  UTILITY 


Figure  7  -  Section  of  Parallel  Conductor  MHD  Magnets 
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and  flux  density  given,  or  minimum  area  and  hence  weight  of  conductor  for  power  or 
j  and  flux  density  given  in  superconducting  or  normal  magnets.  Results  obtained  from 
the  formulas  above  for  cases  of  interest  are  given  in  Table  A-l.  Figure  8  shows  results 
for  dimensionless  power  of  these  cases. 


I 

I 

I 

J 

1 

I 
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Table  A-l 


A/a2 


Figure  8  -  Minimum  Power  in  Optimized  Air  Core 
Electromagnets  for  MHD 
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APPENDIX  B 


MAGNETIC  ENERGY 


The  vector  potential  for  parallel  infinite  rectangular  conductors  as  shown  in  Figure  7, 
Case  4,  due  to  the  sum  of  all  current  elements  j  dx1  dy1  is: 

pi,  f 

Y  '  *  ’ 


(B-l) 


Ay  =  A  =0,  Az=  j  In  [(x-x')2  +  (y-y')2]  dx'dy' 
conductor  area 


4>(x-aa,  y-pa)-<j»(x-a,  y-pa)-<|>(x-aa,  y+pa)+^(x-a,  y+pa) 

+$(xf  aa,  y-pa)-$(xt-a,  y-pa)-$(xf  aa,  y+|3a)+$(>d-a,  y+pa) 

l  /  \  1/2,  2.  2  -1u  2  -1u 

where  ^(u,v)  =  uv  In  (u  +v  )  +  v  tan  —  +  u  ctn  — 


Magnetic  energy  per  unit  length  has  been  found  using  (B-l)  in  the  integration: 

Pa  aa 


r  r  Pa  aa 

I  =  ’  jJj  .J- A  dVol  =  2i  i  I' 

|  2|  coirvoTume  0  J0  -'a 

.2  4 


Az(x,  y)  d  x  dy 


I=2p 1° 

|  JiT 


,a-l.4 ,  ,a-l  2  1  4  2  ,crtK4  .  /CH-1? 

-  (-5—)  ln(-j-)  +2  a  Ina  -  (y-)  In  (-y^ 


-p4  In  p2  -■!  (l-6p2+p^)  In  (1+p2) 


(^l)4  -6^2  pV 


-  i  (a4  -  6  aV+p4)  In  (a^p2) 


In 


(^4  -  6 

2 


/CH-1, 


2.4 


In 


(y4  +P2 


"  (y42+  P2 


8(2^1)  P  (Jl)  -P2  tan'V^) 


+  4  ap  (a^p2)  tan  ^  |3/a 

-8(^1)  p  P2 


fan”1  fi/(y4 


+  4p(l-p2)tan‘1p-4iT(fjl)  p3 
This  result  agrees  with  that  by  Dwight,  Reference  7. 


(B— 2) 
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APPENDIX  C 


PROPAGATION  OF  A  NORMAL  FRONT  IN  A  SUPERCONDUCTING  SYSTEM 


Steady  state  thermal  propagation  of  a  normal  front  along  a  superconductor  may  be  obtained 
by  considering  an  infinite  wire,  cooled  by  a  liquid  at  temperature,  T  ,  with  an  overall 
heat  transfer  coefficient,  h.  The  region  x>  o  is  in  the  normal  state.  Fora  coordinate  sys¬ 
tem  moving  with  the  interface  velocity,  v,  the  heat  transfer  problem  is: 


dT 


dT 


•2P 

i 


ar(kA  3x}  “pcvAir " hs  (Tn  “  V  +  hr ~ °'  x>0 


.  dT  dT 

3T  (kA^f}  "  cvA17  “  hs  ( W  =  °'  x  <  0 

Ts(°)=  Tn  (0)  =  Tc;  Ts(-  oo)  =  To;  ^  (co)  =  0;  k  J-(0)  =  (0) 


>  Eq.  (C-l) 


1  - 


2hsA(Tc-To) 


TT~ 

i  Pe 


h  -  hsA(Tc-To) 

f  ~5 - 

1  P» 


Eq.  (C-2) 


This  result  is  the  same  as  that  given  by  Broom  and  Rhoderick,  Reference  81  When  cooling 
is  not  present  at  the  wire  surface  (h  =  o),  the  bracket  in  Equation  (C-2)  reduces  to  unity. 
Using  Equation  (C-2)  and  average  thermal  properties  obtained  by  combining  winding 
layers,  insulation,  and  secondary  material,  the  smallest  propagation  veloc ity,  that  perpen¬ 
dicular  to  the  layers,  is  calculated  to  be  375  ciry^sec  for  the  magnet  described  here.  Much 
higher  velocities  along  the  wire,  as  might  be  predicted  by  Equation  (C-2),  cannot  be 
obtained  because  diffusion  to  helium,  causing  development  of  film  boiling,  or  in  coils, 
diffusion  to  high  specific  heat  plastic  insulation  increases  effective  pc  in  Equation  (C-2). 


Velocity  along  layers  would  be  somewhat  greater  than  across  layers  because  of  higher 

2 

effective  k.  fhe  result  shows  that  for  2hsA  (T^  -  Tq)>  i  p  ,  propagation  will  not 
occur. 


Riemersma,  Reference  5,  experimentally  investigated  propagation  in  small  solenoids. 

The  results  in  Figure  9  show  thermal  propagation  in  the  A  coil  until  f  equal  0.  25,  when 
external  cooling  and  complete  energy  dissipation  stopped  propagation.  A  similar  coil,  B, 
with  turn  to  turn  shunting  propagated  normalcy  fifteen  times  faster  by  the  coupling  pro¬ 
cess  described  earlier.  The  coil  B  had  1207  turns  of  0.010  inch  Nb-Zr  wire  wound  into  a 

-8 

0. 19  inch  I.  D. ,  0. 5  inch  O.  D. ,  1  inch  long  coil.  The  per  turn  was  2. 13  x  10 

henries;  R.  per  turn,  0.113  ohms;  k,  0.727  turn  to  turn,  and  average  i  (through  the  coil 
1  c 

and  during  the  transient)  was  assumed  10%  higher  than  i^.  We  then  predict  the  time  for 

j  • 

the  coil  to  become  normal,  c 


o  L1  ,l-l£v 

-  Trj-  "  TT')' 


to  be  15  microseconds.  This  result  is 


plotted  in  Figure  9,  whege  it  may  be  seen  that  there  is  good  agreement,  especially  in  the 
rate  when  30%  to  70%  of  the  coil  is  normal.  A  more  accurate  prediction  could  be  obtained 
using  critical  current,  field  and  coil  current  decay  relations  to  eliminate  the  approximation 
of  constant  i^  -  i ^i  in  the  theory.  Also,  improved  correlation  would  be  obtained  by  con¬ 
sidering  slower  transients  in  all  neighboring  turns  and  layers.  The  comparative  experiments 
certainly  demonstrate  that  induction  propagation  exists  and  that  rates  of  the  order  of 
25,000  crr\/sec  may  be  predicted  with  reasonable  accuracy. 
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COIL  a:  1380  TURNS,  670  GAUSS/AMP,  EPOXY  INSULATION. 

COIL  B!  1207  TURNS,  585  GAUSS/AMP ,  SHUNTED  BY 
CONDUCTIVE  INSULATION  (0.01  OHM  CM.) 

COILS  ARE  0.5  IN.  O.D.  x  0.19  IN.  I.D.  x  1.0  IN.  LONG,  0.010  IN. 
WIRE,  0.0005  IN.  THICK  INSULATION. 


Figure  9  -  Decay  of  Magnet  Field  and  increase  of  Resistance 

Upon  Entry  to  Normal  State  of  Shunted  and  Unshunted  Test  Coils 
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PERCENT  NORMAL  RESISTANCE 


APPENDIX  D 


MAGNETIC  STRESS  AND  FORCE 

Stress  resulting  from  j_  x  B  body  force  per  unit  volume  is  easily  found  considering  that  B 
=  \7x  A  and  that  in  Figure  7,  Case  4,  boundaries  x  =  a  and  y  =  pa  have  zero  stress. 


>/y)  =  -i QJ  By(x,y)dx  =  j£ 


Pa 


%(x'y)  =  -i  0f  Bx(x/y)dy  =  i 

J  y 


A  (x,  y)  -  A  (a,  y) 


A  (x,  y)-A  (x,  pa) 


l  Eq.  (D— 1 ) 


Net  force  F  separating  conductors  may  be  found  either  by  integrating  compressive  stress 
at  the  outer  edge  or  using  the  principle  that  the  variation  in  total  magnetic  energy,  T, 
during  a  relative  displacement  of  conductors  at  constant  current,  is  equal  the  work  by 
mechanical  forces  done  by  the  current  system.  See  for  example  Dwight,  Reference  9. 
Keeping  coil  sections  constant,  especially  (a-l)a. 


2  2 
a  +3 
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«^>  +  B2  J 
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T~J 


+  3p2  in 
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<¥>  *P2 
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2  . 
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-  In  |l+p2j  +  ^P  (3a^  “  P^)  ton  ^  P^a 

2 

+  (3-p2)  tan  P  -  2  (3  -  p2)  tan  ^  p/(— *-) 


Similarly  end  force  is  found  as: 

_  1  dT  _ 1  T 

e  ~T  3T  T 


Eq.  (D-2) 


Eq.  (D-3) 
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APPENDIX  E 


VAPOR  COOLING  OF  THERMAL  CONDUCTION 
PATHS  AND  ELECTRICAL  LEADS 


Heat  transfer  in  a  thermal  and  electrical  conductor  cooled  by  a  mass  flow  m  from  a  boiling 
fluid  at  one  end  of  the  conductor  is  summarized  below: 

.2 


-  MT  -  Tv)  -  0 


T  (0)  =  T  (0)  =  T  ,  T  (x  )  =  T 
v  o  o  a 


For  constant  properties: 


T  "  T~  “  (Ta  '  V  |  <Xo}  "  *1  (x) 


G+?\  I  eEx-l 

ft  hr 
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/v  I  W  1  I  U 

*2(x)_  w  e  '  T+r 


G-El  /  e'Fx-l 
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G+Fl  Ex  IG-E 


F  _  hs  .  (2m  Co)2  , 

E  r +  rax  1 


,  F  = 


-Fx  . 
e  -1, 
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2fal 


1  + 


(2mc  ) 


hskA 


>  Eq.  (E-l) 


)  Eq.  (E-2) 


me 


f+1  hG=l^ 


The  heat  which  is  conducted  t^  boiling  fluid  is 


1  P„ 


„  ♦,<*«> 

Q°  ■"♦jTg - 


Eq.  (E-3) 


When  h  is  considered  infinite,  T  =  T  ,  the  result  becomes  from  Equation  (E-3): 
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Eq.  (E-4) 
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1  pe  x  2 

me  (T  -T  )  — r-2  o  V  p  K 


Q  =-P-q— °- 


mc  x 


■+  -r- 


mc 


ta 


If  Qq  is  the  only  heat  flux  vaporizing  liquid  then  it  is  also  mh^.  Then  where  |  is  zero: 
.  kA 


c(T  -T  ) 

m  = -  In  (1  +  ■ P I — ■ — ) 

c  x  h,  ' 

p  o  fg 


Eq.  (E-5) 


M 


inimizing  Qq  in  Equation  (E-4)  with  respect  to  x JA  for  an  electrical  lead,  consider¬ 


ing  me  x  /kA  larger  than  one  gives: 
P  ° 


x 

o 

~K 


optimum  lead  me 


mcn(T  ) 

p  a  o 

“2 - 

1  P 

'  A 


Eq.  (E— 6) 


The  minimum  value  of  mc^x ^kA  results  when  m  is  just  Q^/h^.  And  using  Equations 
C  (Ta  -V 

(E-4)  and  (E-6)  is  1  +  Pi  ° - Py  which  is  considerably  greater  than  one  for  helium  to 

fg 

nitrogen  or  room  temperature.  Thus  the  optimization  is  rather  meaningless  since  the  heat 
2 

flux  Qq  is  practically  i  k/ me  for  a  range  of  values  of  xq/A,  if  m  results  independently 

from  a  number  heat  flux  sources  other  than  leads. 
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VAPOR 

FLOW  WARM 


T  w 

0  COLD 


Figure  10  -  Vapor  Cooling  of  Conduction  Paths  and  Electrical  Leads 
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APPENDIX  F 


VAPOR  COOL-DOWN 


Cooling  a  mass  M  with  boiling  helium  from  a  temperature  to  requires  a  mass  m  of 
helium: 

T, 


=  M 

mboiling  TTT 
f9 


/ 


cdT 


Eq.  (F-l) 


But  if  vapor  is  also  used  to  cool  an  isothermal  mass  by  passing  through  an  isothermal  heat 
exchanger  at  temperature  T,  a  smaller  amount  of  helium  is  needed. 

T, 

ca  i 

Eq.  (F-2) 


m  ,=  M 

vapor  cool 
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/  V+(W°) 


hA 


T  »-T 

vapor  out  o  _  e  _  | 
^  To 


wall 


me 

■e  p 


Eq.  (F-3) 


For  helium,  h^  -  eCp  T^,  is  approximately  zero.  Most  metals  have  specific  heats  which 
vary  as  c  =  c  (T/T  )n,  2<n<3. 


o  '  o 

m 


vapor 


cool  _  %  |  n+1 


m. 


boiling 


ec  T.  In 
P  1 


Eq.  (F-4) 


For  heat  exchange  efficiency  of  0.75,  cooling  metal  from  nitrogen  temperature,  77  K, 
to  4. 2°K  produces  a  ratio  of  9  to  11  from  (F-4)  while  a  factor  of  7  is  obtained  for 
plastics. 


Consider  that  the  body  M  is  not  isothermal  but  made  up  of  differential  masses  M  dx/ x^, 
that  specific  heat,  c,  is  constant,  that  helium  vapor  flow  starts  at  x  =o,  that  there  is 
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always  equilibrium  with  the  vapor  and  the  masses  do  not  conduct  heat  to  each  other. 
Then: 


me  x 
P  o 
Me 


§I+fI  '  0.  T  (0,  t)  =  Tq,T(x,o)  *  T, 


-1 
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T  \  — 
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_  .  c  x 

(T,  T  ) - —  p 

i-o  me  x 

e  p  o 


Eq.  (F-5) 
Eq.  (F-6) 


,  P  P 

This  result  shows  that  mass  and  vapor  temperatures  are  a  step  function  up  from  T 


me 


to  Tj  at  ~  Thus  all  enthalpy  of  the  vapor  between  TQand  T^  is  used  until 

the  last  mass  at  xq  is  cooled.  A  maximum  possible  reduction  of  m  from  that  by  Equation 
(F-l)  is  therefore  1  +  c^  (T^  -  TQ)/h^,  or  for  helium  to  nitrogen  temperature,  19.6. 
Isolating  system  parts  and  cooling  down  with  a  series  path  for  vapor  has  some  advantage. 
But  isothermal  mass  cool  down  with  vapor  and  vapor-mass  equilibrium  (e  =  1)  gives  factors 
only  25%  to  33%  less  because  materials  have  large  specific  heats  at  high  temperatures. 
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APPENDIX  G 


LEVEL  GAGES 


Convection  heat  transfer  coefficients  from  a  resistance  element  to  vapor  are  much  less 

than  in  boiling  in  liquid.  Temperature  rise  causes  a  proportional  resistivity  change  in 

vapor  and  a  variation  of  resistance  with  level.  A  conductor  of  area  A,  perimeters  s 

for  convection  and  having  temperature  coefficient  of  resistivity  a  and  effective  heat 

conduction  area  A1  when  there  is  support  structure,  immersed  in  liquid  at  T  from  x 

o  o 

to  |  ,  but  vapor  at  Tq  up  to  xq  is  described  by: 
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Eq.  (G-l) 
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Eq.  (G-3) 


Gage  sensitivity  a6  is  also  a  (gage  power)/hS.  Therefore  it  is  important  to  have  high  and 
small  wire  perimeter  to  minimize  power.  Ten  feet  of  0.00124  inch  diameter  copper  wire 
wound  on  a  mylar  ribbed  coil  form  serves  as  a  nitrogen  gage  consuming  0.01  watts.  The 
helium  gage  is  a  string  of  eight  carbon,  1/10  watt,  27  ohm  resistors  consuming  less  than 
0. 001  watts.  Sharpness  parameter  8,  is  adequate  for  usefulness  for  both  gages  that  are 
shown  In  Figure  11,  although  there  is  not  linearity.  The  measuring  circuit  used  is  similar 
to  that  used  by  Maimoni,  Reference  10,  who  gives  a  good  discussion  of  resistance  level  gages. 
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Helium 


Nitrogen 


Figure  11  -  Helium  and  Nitrogen  Level  Gages 
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NOMENCLATURE 


half  width  between  magnet  conductors 

section  area  of  magnet  winding,  thermal  path,  etc. 

n'th  Bernoulli  number 

magnetic  flux  density 

specific  heat  of  magnet,  etc. 

specific  heat  of  vapor 

heat  exchanger  effectiveness 

fraction  of  coil  that  is  normal  (resistive) 

force 

convection  heat  transfer  coefficient 
latent  heat  of  vaporization 
magnetic  field  intensity 
current 

current  density 

thermal  conductivity  or  magnetic  coupling  coefficient 
MHD  duct,  magnet  or  gage  wire  length 
inductance 

inverse  Laplace  transform  operator 
mass  of  coolant 

mutual  inductance  or  mass  to  be  cooled 
number  of  shunted  coil  sections 
number  of  layers  in  coil 
Laplace  transform  variable  • 

generated  power 
electrical  power 
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Q  = 

r  = 


R 


s 


T 

u  = 

V  = 

V 

n 

W 

x,  y,  z  = 

a  = 

a  = 
r 

°Va2  = 
(3  = 

y 

e  = 


*n 

^e 

X 


H 

p« 

P 


a 


a 

x 

T 


/ 


a 

y 


heat  removed  at  temperature  T 

radius 

resistance 

gas  constant  of  helium 
perimeter 

absolute  temperature 

plasma  velocity 

velocity  of  normal  front 

voltage  across  normal  fraction  of  coil 

system  weight  to  power  ratio 

coordinates 

dimensionless  distance  to  outside  of  magnet 
temperature  coefficient  of  resistance 
accommodation  coefficients 
dimensionless  half  height  of  magnet 
ratio  of  specific  heats 
emissivity 

ratio  of  coefficient  of  performance  to  Carnot  cycle  COP 

electrical  efficiency 

ratio  conductor  volume  to  total  volume 

magnetic  permeability  of  free  space 

electrical  resistivity 

density 

plasma  conductivity  or  Stefan-Boltzmann  constant 

stress 

time 
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